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ABSTRACT 

Wc present a study of Spitzer /IRAC and X-ray active galactic nuclei (AGNs) selection techniques 
in order to quantify the overlap, uniqueness, contamination, and completeness of each. We investigate 
how the overlap and possible contamination of the samples depends on the depth of both the IR 
and X-ray data. We use Spitzer /IRAC imaging, Chandra and XMM-Newton X-ray imaging, and 
spectroscopic redshifts from the PRism MUlti-object Survey (PRIMUS) to construct galaxy and 
AGN samples at 0.2 < z < 1.2 over 8 deg 2 . We construct samples over a wide range of IRAC 
flux limits (SWIRE to GOODS depth) and X-ray flux limits (10 ks to 2 Ms). We compare IR- 
AGN samples defined using the IRAC color selection of Stern et al. (2005) and Donley et al. (2012) 
with X-ray detected AGN samples. For roughly similar depth IR and X-ray surveys, we find that 
^75% of IR-selected AGN are also identified as X-ray AGN. This fraction increases to ^90% when 
comparing against the deepest X-ray data, indicating that only ~10% of IR-selected AGN may be 
heavily obscured. The IR-AGN selection proposed by Stern et al. suffers from contamination by 
star-forming galaxies at various redshifts when using deeper IR data, though the selection technique 
works well for shallow IR data. While similar overall, the IR-AGN samples preferentially contain more 
luminous AGN, while the X-ray AGN samples preferentially contain lower specific accretion rate AGN, 
where the host galaxy light dominates at IR wavelengths. The host galaxy populations of the IR and 
X-ray AGN samples have similar restframe colors and stellar masses; both selections identify AGN in 
blue, star-forming and red, quiescent galaxies. 

Subject headings: galaxies: active — galaxies: evolution — infrared: galaxies — X-rays: galaxies 



1. INTRODUCTION 

Understanding the nature and role of active galactic 
nuclei (AGNs) is crucial for understanding both the ac- 
cretion history of the Universe as well as galaxy evolution. 
There is mounting observational evidence that there is a 
connection between black hole growth and galaxy growth. 
This is shown both by the tight correlation between black 
hole mass and galaxy bulge mass (e.g. Magorrian et al. 
1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000) 
as well as the similar evolutionary history of star for- 
mation and AGN activity through cosmic time (e.g. 
Boyle & Terlevich 1998; Silverman et al. 2008; Aird et al. 
2010). In order to understand and characterize AGNs, 
one must be able to identify a complete AGN sample, 
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with full knowledge of any underlying biases or contam- 
ination in the sample. One can then better determine 
which processes (e.g. secular evolution, mergers, envi- 
ronment) are the dominate fueling mechanisms in the 
growth and evolution of AGNs. 

Deep X-ray surveys provide a reliable means of select- 
ing AGN, in that they introduce few false positives, as 
the AGN light generally outshines light from even highly 
active star-forming galaxies at X-ray wavelengths. How- 
ever, high-column density gas (Njj > 10 23 cm -2 ) absorbs 
X-rays, such that X-ray surveys may fail to identify the 
most heavily absorbed AGNs. While the exact fraction of 
obscured AGNs is not yet fully known, obscured AGNs 
likely represent a large fraction of the total AGN pop- 
ulation at all luminosities. For example, Treister et al. 
(2004) predict that 25% to 50% of AGNs are obscured, 
even at high luminosities (Lx > 10 44 ergs _1 )(see also 
Ballantyne et al. 2006), while Gilli et al. (2007) uses the 
cosmic X-ray background to predict that Compton-thick 
AGNs are four times as numerous as unobscured AGNs 
at low luminosity (Lx < 10 43 ' 5 ergs -1 ). Only with ex- 
tremely deep X-ray data (>1 Ms) docs it become pos- 
sible to detect and characterize the heavily obscured, 
moderate luminosity AGNs (e.g. Georgantopoulos et al. 
2009; Brightman & Ueda 2012) that would otherwise be 
missed in shallower X-ray surveys. However, even the 
deepest X-ray surveys may still fail to identify low lumi- 
nosity AGNs with moderate to heavy obscuration. Ad- 
ditionally, current deep X-ray surveys cover at most 
~1 deg 2 of sky, which limits the size of the resulting 
AGN samples. Ideally, both wide and deep X-ray sur- 
veys are required to identify statistically large, relatively 
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complete AGN samples. 

The use of mid-infrared (MIR) emission to identify 
AGNs began in the 1970s with the advent of sensi- 
tive IR detectors (Low & Klcinmann 1968) and contin- 
ued with the Infrared Astronomical Satellite (IRAS) and 
the Infrared Array Camera (IRAC Fazio et al. 2004) on 
board Spitzer. High-energy radiation from the AGN is 
reprocessed by dust near the AGN and re-radiated at 
MIR wavelengths. Luminous AGNs display a red MIR 
power- law SED, which is dominated by thermal emis- 
sion from hot dust (Neugebauer et al. 1979; Elvis et al. 
1994; Rieke & Lebofsky 1981). Emission at MIR wave- 
lengths can also potentially be used to detect Compton- 
thick AGNs that may be missed by deep X-ray surveys 
(e.g. Ivison et al. 2004; Lacy et al. 2004; Stern et al. 2005; 
Alonso-Herrero et al. 2006; Polletta et al. 2006). It has 
also been found that objects with similar MIR luminosi- 
ties can have vastly different radio, optical, UV and soft 
X-ray luminosities (Mushotzky 2004), such that heav- 
ily obscured AGNs can be identified in the MIR as it 
is relatively insensitive to obscuration. Therefore IR- 
AGN identification is a potentially powerful tool that is 
sensitive to both obscured and unobscured AGNs, with- 
out requiring time-consuming deep X-ray data. Build- 
ing on the above findings, a variety of selection tech- 
niques have been developed to identify large samples of 
AGNs from MIR imaging data. Within the MIR wave- 
band there arc many techniques to select infrarcd-AGN 
(IR- AGN) , each taking advantage of the unique colors of 
AGNs. Stern et al. (2005) and Lacy et al. (2004) pro- 
pose IRAC color-color cuts that were designed using 
shallow IRAC surveys and which effectively select lu- 
minous AGNs at low redshift. These selections identify 
the AGN population and separate them from the much 
larger galaxy population based on their characteristic 
MIR properties. Alternatively, the IR SED can be fit 
using multiple bands and the associated photometric er- 
rors to estimate the probability that the source has a fea- 
tureless power-law continuum (e.g. Alonso-Herrero et al. 
2006; Polletta et al. 2006; Donley et al. 2007). Such tech- 
niques identify AGNs within similar regions in color-color 
space but are generally more reliable and lead to smaller 
samples. 

Recent work using deeper IRAC and X-ray surveys (e.g. 
Barmby et al. 2006; Donley et al. 2007; Cardamone et al. 
2008; Park et al. 2010) has begun to investigate these 
IR-AGN selections beyond the shallow surveys for which 
they were designed. Barmby et al. (2006) use deep IR 
and X-ray data in the Extended Groth Strip (EGS) to 
find that X-ray AGNs have a wide range of MIR col- 
ors, suggesting that there is no single method that will 
be able to identify a complete AGN sample. Likewise, 
Park et al. (2010) find that a majority (78%) of the X- 
ray AGN sources are not detected by the power-law AGN 
selection method. Together these studies suggest that 
there is a population of AGNs missed by IR-AGN se- 
lection techniques even at shallow X-ray and IR survey 
depths. Recently, Donley et al. (2012) revisit the Stern 
et al. and Lacy IRAC color-color selections, which were 
defined with relatively shallow surveys, and the power- 
law selection, which depends on both the photometry 
and (generally underestimated) photometric errors, and 
define a new IRAC color-color selection criteria, which 
they claim is more reliable than the previous IR selec- 



tion criteria. While each of these IR-AGN criteria select 
some fraction of the total underlying AGN population, 
Barmby et al. (2006) suggest that no proposed IR-AGN 
color selection will identify all AGNs, given the wide 
range of spectral shapes exhibited by X-ray sources in 
the EGS. 

With some possible issues, using the MIR to se- 
lect AGNs is still an effective way to identify heavily 
obscured/Compton-thick AGNs that are an important 
population missed by X-ray surveys. In order to build 
more complete AGN samples it has become common 
to use multiple waveband selection techniques in the 
IR through X-ray. For example, Hickox et al. (2007), 
and Assef et al. (2010) combine different waveband se- 
lection techniques to probe the properties of the underly- 
ing AGN population. Combining both X-ray AGN sam- 
ples and IR-AGN samples should provide a better census 
of the underlying AGN population. However, there has 
been relatively little study of the overlap and uniqueness 
of each AGN selection technique. 

Comparing samples selected using different techniques 
could provide insight into the AGNs that are missed us- 
ing different selection techniques. Similarly, it is impor- 
tant to understand the properties of AGNs and their host 
galaxies selected in each way. 

Here we aim to quantify the overlap and uniqueness of 
various IR-AGN selection techniques compared to X-ray 
AGN selection. We use X-ray and IR data of varying 
depths in the CDFS, COSMOS, ESI and XMM fields 
covered by the PRIMUS redshift survey. We use these 
datasets to investigate how the selected AGN population 
depends on the selection technique and the depth of the 
data, and we investigate the completeness, contamina- 
tion, and uniqueness of X-ray versus IR-AGN selection. 
The paper is organized as follows. In §2, we present the 
relevant multi- wavelength datasets. In §3, we detail the 
different AGN selection techniques. In §4, we compare 
the number densities and overlap between X-ray AGN 
and IR-AGN selection techniques, as a function of sur- 
vey depth. In §5, we investigate contamination of IR- 
AGN samples. In §6, we compare the AGN and host 
galaxy rest-frame properties of X-ray and IR selection 
techniques. We discuss our results in §7 and conclude 
in §8. Throughout the paper we assume a standard flat 
ACDM model with fi m = 0.3, fl A = 0.7, and H = 72 km 
s _1 Mpc" 1 . 

2. DATA 

For this study, we use multi-wavelength data from 
the COSMOS, ELAIS-S1 (ESI), XMM-LSS(XMM) and 
Chandra Deep Field South (CDFS) fields. All of these 
fields have optical photometry, Spitzer IR imaging, spec- 
troscopic redshifts from PRIMUS, and X-ray data from 
Chandra and XMM-Newton. We describe these datasets 
in detail in Sections 2.1, 2.2, and 2.3 below. Figure 1 
shows the IRAC (orange boxes), PRIMUS (blue), and 
X-ray coverage (green) for each of the fields. 

2.1. Spitzer IR Data 

In the CDFS, ESI, and XMM fields, we use ex- 
isting Spitzer Infrared Array Camera (IRAC) imaging 
from Data Release 2 (DR2) of the Spitzer Wide- area 
Infrared Extragalactic Survey (SWIRE) (Lonsdale et al. 
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Figure 1. Map of the IR (orange box), PRIMUS (blue), and 
X-ray (green) coverage in the CDFS, COSMOS, ESI, and XMM 
fields. Ten thousand random objects detected in all four IRAC 
bands are shown with grey points. Within CDFS the red outline 
shows the SIMPLE region, which has deeper IR data, and the solid 
red polygon shows the GOODS region, which has the deepest IR 
and X-ray data. The SCOSMOS footprint is outlined in red in 
the COSMOS field, and the deeper X-ray footprint in the central 
0.9 deg 2 of the COSMOS field is outlined in black. 



2003). Within the shallow CDFS-SWIRE (orange outline 
in Figure 1), we include deeper source catalogs within 
the CDFS-SIMPLE region (solid red outline) and even 
deeper data in CDFS-GOODS region (solid red poly- 
gon), both taken from version 3.0 of the Spitzer IRAC 
/ MUSYC Public Legacy in E-CDFS (SIMPLE) survey 
data release (Damen et al. 2011). In the COSMOS field 
we use existing IRAC imaging from the S-COSMOS Sur- 
vey (Sanders et al. 2007). 

In fields with SWIRE coverage, we adopt source cat- 
alogs provided by the SWIRE team as part of their 
SWIRE Data Release 2 12 . These catalogs were gener- 
ated by performing source detection in each of the four 
IRAC channels and merging the source lists (see below 
for more details). For the CDFS-SIMPLE region with 
deeper IRAC data, we compared the measured fluxes for 
sources in both the SWIRE catalog and the SIMPLE cat- 
alog (Damen et al. 2011). Based on this comparison, we 
find that the SIMPLE fluxes need to be scaled by 2.02, 
1.96, 1.69, and 1.56, in IRAC channels 1, 2, 3, and 4 re- 
spectively to match the measured fluxes from the CDFS- 
SWIRE catalog. This multiplication factor accounts for 
the total flux correction applied by the SWIRE team com- 
pared to the 4.0" aperture used by the SIMPLE team and 
a small zero point difference. The errors are also multi- 
plied by the above factors, and agree with the SWIRE er- 
rors after accounting for the difference in exposure times. 
The SIMPLE catalog includes our deepest IRAC data in 
the CDFS-GOODS region. 

For the S-COSMOS IRAC data, we reproduce the 
SWIRE source detection procedure (as outlined in the 
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Figure 2. Number density (number of sources per logarithmic 
flux bin per square degree) detected in all four IRAC channels as 
a function of IRAC 5.8 fim flux. We show number densities for 
our three survey depths: SWIRE (orange squares), COSMOS (purple 
pluses), and GOODS (light blue crosses). The vertical colored lines 
indicate our chosen 5.8 fim flux limit for each survey depth, which 
is set at a flux brighter than the observed turn over in the number 
counts (see Table 1). The solid black line shows the combined 
number density of sources across all fields that reach our fixed 
5.8 fim flux limits. We also show the Stern et al. , Donley et al. , 
and power-law IR-AGN total number densities from the combined 
sample (grey lines). 

SWIRE DR2 documentation) to ensure we measure ro- 
bust fluxes and errors using a consistent technique. We 
downloaded the mosaic images in each IRAC channel 
from the NASA/IPAC Infrared Science Archive 13 (May 
2007 release). We use Astronomical Point source Ex- 
tractor (APEX) pipeline in MOPEX to perform source 
detection in the mosaic for each channel on an individ- 
ual basis. We combined the individual catalogs using 
the BANDMERGE procedure in MOPEX. The merged 
sources in the final catalog are required to be detected 
with Signal-to-Noise Ratio (SNR) > 10 in channel 1 and 
SNR > 5 in channel 2. We follow the SWIRE handbook 
and apply a "Coverage-SNR" flux limit, which is a factor 
proportional to the inverse of the square root of the cov- 
erage (C), the number of times that the object's location 
was scanned by IRAC in a given channel. This threshold 
limits flux measurements to where the flux of a detected 
object is large compared to the statistical fluctuations 
in the background sky in a given channel. Due to the 
deeper coverage of Spitzer IRAC scans in the COSMOS 
field, we require that the coverage threshold to be 3, 3, 
16, and 50 /iJy x \/4/C for channels 1, 2, 3, or 4, respec- 
tively to match the SWIRE fields. For the majority of 
sources, our flux measurements are similar to those in the 
S-COSMOS public catalog, although the public catalog 
tends to have larger errors for all objects. 
We combine our various fields and datascts into three 
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"IR surveys depths" with data of comparable depth: 

1. GOODS, the deepest area of IRAC imaging within 
the CDFS held; 

2. COSMOS, including the S-COSMOS data and the 
SIMPLE coverage in the CDFS field that reaches 
comparable flux limits; 

3. SWIRE, which includes the entire area of the XMM 
and ESI fields as well as the large area in the CDFS 
field with SWIRE coverage. 

In Figure 2 we show the number density of sources for 
each survey depth as a function of 5.8 fj,m flux. We re- 
quire a source to be detected in all four IRAC bands, 
which is necessary for our analysis in this paper. We 
set a 5.8 fim flux limit for each IR survey depth that is 
just above the point where the number counts start to 
turn over (vertical colored lines), indicating incomplete- 
ness in the IRAC detected samples. These flux limits are 
given in Table 1. Our SWIRE IR survey depth includes 
the entire XMM, ESI and CDFS fields (21.55 deg 2 in to- 
tal); our COSMOS IR survey depth includes COSMOS and 
the SIMPLE region of the CDFS field (3.05 deg 2 ). Our 
deepest IR survey depth (GOODS) includes only the small 
GOODS region of the CDFS field (0.06 deg 2 ). In the rest 
of this paper we adopt these limits and the correspond- 
ing area coverages when analyzing data to a given IR 
survey depth. The solid black line in Figure 2 shows the 
overall number density combining all of our fields that 
probe above the flux limit for each IR survey depth. 

2.2. PRIMUS Data 

PRIMUS (Coil et al. 2011) is the largest intermediate- 
redshift, faint-galaxy spectroscopic survey performed to 
date and covers ~9 deg 2 in seven different fields with 
existing deep multi-wavelength imaging. Wc obtain low- 
resolution (R ~ 40) spectra for ^300,000 objects from 
the IMACS instrument on the Magellan I Baade 6.5m 
telescope, targeting 80% of galaxies with i < 22, with a 
statistically complete sample to a depth of i ~ 23. By 
fitting galaxy, broad-line AGN, and stellar spectral tem- 
plates to low-resolution spectra and optical ground based 
photometry we have measured ^120, 000 robust redshifts 
with a precision of cr 2 /(l + z) < 0.5% and catastrophic 
outlier rate of < 3% (Az/(1 + z) > 0.03). We classify 
objects as galaxies, broad-line AGN or stars based on 
the x 2 °f the best, template fits. We derive K-corrections 
(Blanton & Rowcis 2007) from the photometry. For fur- 
ther details of the survey design, targeting, and data see 
Coil et al. (2011); for details of the data reduction, red- 
shift confidence and precision, and completeness see Cool 
et al. (in preparation). 

Here we use PRIMUS redshifts between 0.2 < z < 1.2 
with high confidence quality flags (Q > 3, see Coil et al. 
(2011).) In the CDFS field, the PRIMUS science ob- 
servations cover part of CDFS-SWIRE but are disjoint 
from the CDFS-SIMPLE and CDFS-GOODS area. We 
therefore also include the PRIMUS CDFS calibration 
field, which overlaps with the CDFS-GOODS and CDFS- 
SIMPLE areas, to increase the number of deep IRAC and 
X-ray detected objects with PRIMUS redshifts for Sec- 
tion 6. We thus have a total of 31,998 galaxies and broad- 
line AGNs with high-confidence PRIMUS redshifts in 



the COSMOS, CDFS, ESI, and XMM fields. Figure 1 
shows the overlaps between the PRIMUS footprint and 
the IRAC and X-ray coverage. Table 2 gives the areas of 
the total and overlapping regions of each data set. 

2.3. X-ray Data 

We have compiled published X-ray source catalogs 
based on existing Chandra and XMM-Newton X-ray 
data in the COSMOS, CDFS, ESI, and XMM fields 
(See Aird et al. (2012) for details). In all of the 
fields we use the likelihood ratio matching technique 
(e.g. Sutherland & Saunders 1992; Ciliegi et al. 2003; 
Brusa et al. 2007; Laird et al. 2009) to identify counter- 
parts for each of the X-ray sources in the PRIMUS tar- 
geting optical photometry. Wc also apply the likelihood 
ratio technique to assign secure IRAC counterparts to 
the X-ray sources in the entire area with both X-ray and 
Spitzer IRAC coverage. Candidate counterparts are iden- 
tified within 5", although over 80% of likelihood ratio 
matched IRAC counterparts are within 2" of their X-ray 
counterpart. The likelihood ratio matching technique ac- 
counts for optical, IR, and X-ray positional uncertainties, 
the probability of having a counterpart with a given mag- 
nitude, and the probability of a spurious match. Where 
multiple counterparts exist, we choose the match with 
the highest likelihood ratio and restrict to "secure" coun- 
terparts with likelihood ratios >0.5. 

Within the CDFS field, we use the Luo et al. (2008) 
2 Ms CDF-S X-ray source catalog, one of the deep- 
est Chandra survey to date, and the four overlapping 
250 ks Chandra point source catalogs of Lehmcr et al. 
(2005). The 2 Ms catalog covers all of the CDFS-GOODS 
IRAC coverage (0.058 deg 2 ), and part of CDFS-SIMPLE 
(0.25 deg 2 ) and CDFS-SWIRE (0.31 deg 2 ) coverage, to 
a depth of /2-8koV ~ 5.5 x 10~ 17 erg s~ 1 cm~ 2 . The 
250 ks catalogs surrounds the CDFS-GOODS field, and 
covers more of the CDFS-SIMPLE (0.068 deg 2 ) and 
CDFS-SWIRE (0.12 deg 2 ) area to a depth of / 2 - 8 kcV ~ 
6.7 x 10" 16 erg s^cm -2 . 

The COSMOS field was observed with XMM-Newton 
(Hasinger et al. 2007) over the entire 2 deg 2 to a depth of 
/2-iokcV ~ 3 x 10~ 15 erg s _1 cm -2 and with much deeper 
Chandra data reaching /2-iokcV ~ 8 x 10~ 16 erg s _1 cm" 2 
for the central ~0.9 deg 2 (Elvis et al. 2009). There are a 
total of 2,769 X-ray sources in the area with both IR 
and X-ray coverage, and 2,075 X-ray sources with X- 
ray, PRIMUS, and IR coverage; 383 of these have robust 
PRIMUS redshifts and classifications. 

We use the Puccetti et al. (2006) point source catalog 
in the ESI field from a mosaic of four partially overlap- 
ping XMM-Newton pointings which reached a depth of 
/2-iokcV ~2x 10~ 15 erg s _1 cm -2 covering 0.52 deg 2 of 
the 0.9 deg 2 of PRIMUS and 2.16 deg 2 with IR coverage. 
Within the PRIMUS and IR areas (see Figure 1) there 
are 371 and 395 X-ray sources, respectively; there are 73 
X-ray sources with secure PRIMUS redshifts. 

The XMM-LSS(XMM) X-ray data are from the 
Subaru/ XMM-Newton Deep Survey (SXDS) (Ueda et al. 
2008), which contains 7 deep XMM-Newton pointings, 
and the XMM-LSS X-ray survey (Pierre et al. 2007) 
which contains 45 pointings. The combined catalog 
contains X-ray sources to a depth of /2-iokcV ~ 2 x 
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and covers 2.16 deg 2 of PRIMUS 



" 15 erg s 'cm 2 , 
and 2.16 deg 2 of the IR area. There are 1,840 and 4,466 
sources in the PRIMUS and IR overlap and 264 ob- 
jects with secure redshifts. 

3. AGN SAMPLE SELECTION 

In this section we define our AGN samples based on two 
IRAC color selection techniques and X-ray detections. 
No single selection technique can identify a complete par- 
ent sample of all AGNs; instead, we must study how 
the overlap between our different samples vary with the 
depths of the observations to determine selection biases 
and incompleteness effects. In Section 3.1, we define our 
Stern et al. IR-AGN sample using the MIR color selec- 
tion criteria from Stern et al. (2005). In Section 3.2, we 
define our Donley et al. IR-AGN sample using the MIR 
color criteria from Donley et al. (2012). In Section 3.3, 
we discuss identification of IR-AGNs by fitting a power- 
law to the MIR photometry (e.g. Alonso-Herrero et al. 
2006; Polletta et al. 2006; Donley et al. 2007) and the 
reasons we choose not to use this technique in the re- 
mainder of this paper. We describe our X-ray selected 
AGN sample in Section 3.4 and explain how we calculate 
completeness weights to account for the varying sensitiv- 
ity of the X-ray observations. 

3.1. Stern Color Selection 

Objects that are detected in all four IRAC bands are 
defined to be Stern et al. selected IR-AGN if they have 
IRAC colors such that they lie within the following region 
in color-color space: 



([5.8] - [8.0]) > 0.6 and 

([3.6] - [4.5]) > 0.2 • ([5.8] - [8.0]) + 0.18 

([3.6] - [4.5]) > 2.5 • ([5.8] - [8.0]) - 3.5 



and 



(1) 
(2) 
(3) 



We apply this selection technique for samples reaching 
our three different "IR survey depths" defined in Section 
2.1 above. We identify a total of 1,137, 10,343, and 6,687 
Stern et al. IR-AGN in our GOODS, COSMOS, and SWIRE 
depth surveys respectively. We show the number density 
of Stern et al. selected objects with X-ray coverage in 
Figure 2. 

3.2. Donley Color Selection 

We also select IR-AGN samples using the IRAC color 
criteria presented by Donley et al. (2012). This color- 
selection technique was designed to limit contamination 
by star-forming galaxies, especially at high redshift, but 
still be both complete and reliable for the identification 
of luminous AGNs. We require that objects are detected 
in all four bands, and have IRAC colors such that they lie 
within the following region in IRAC color-color space: 



x = logio 



/5.8/in 



V = log 



10 
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x>0.08 and y > 0.15 and 
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y<(1.21 x x) +0.27 and 
and / 5 

h.Oum > fb 
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Like in Section 3.1, we apply this selection method to 
our three different IR survey depth samples. Wc identify 
a total of 212, 2,657, and 4,151 Donley et al. IR-AGN in 
our GOODS, COSMOS, and SWIRE depth surveys respectively. 
The number density of Donley et al. selected objects with 
X-ray coverage is shown in Figure 2. 

3.3. Power-law Selection 

We also investigated the identification of IR- 
AGN samples using a power-law selection tech- 
nique Alonso-Herrero et al. (2006); Polletta et al. (2006); 
Donley et al. (2007). This technique identifies IR-AGNs 
by fitting a power-law function to the IRAC fluxes, in 
contrast to the color selection techniques such as the 
Stern et al. and Donley ct al. methods described above. 
Sources dominated by a power-law IR-AGN have a nega- 
tive slope a < —0.5 such that their vF v flux increases 
with wavelength (Alonso-Herrero ct al. 2006). Follow- 
ing Donley et al. (2007), wc identify sources with four 
detected IRAC flux measurements as Power-law AGNs 
when they are well fit by a power law of slope a < —0.5, 
where /„ a c". We deem the power-law fit as good if 
P x > 0.1 (probability that the fit would have a x 2 greater 
than the found \ 2 ). We calculate the P x and fit param- 
eters using the linf it procedure in IDL. 

As in Sections 3.1 and 3.2, we apply the power-law se- 
lection technique to our SWIRE, COSMOS, and GOODS IR 
survey depths. To properly simulate a shallow SWIRE- 
like survey in our COSMOS depth fields (COSMOS and 
CDFS-SIMPLE), wc apply not only the 5.8 fim flux lim- 
its but also scale up the errors in the flux measurements 
to accurately represent the characteristics of the shal- 
lower data. Specifically, we scale the errors in the fluxes 
by the square root of the differences in the exposure times 
for the COSMOS and SWIRE depth data. We verify that the 
median of our scaled errors matches the median error in 
the SWIRE depth fields. 

Figure 2 shows the number density of IR-AGNs identi- 
fied using the power-law selection technique as a function 
of 5.8 [im flux. Unlike the Stern et al. and Donley et 
al. samples, the number density of the Power-law sample 
is not a smooth function of IR flux. Donley ct al. (2012) 
shows that this is mainly due to the flux errors associated 
with each survey. While we have ensured similar data re- 
duction and consistent error estimates for the different 
fields and samples, the uncertainty in a flux measure- 
ment will depend on the depth of the data. At a given 
flux, deeper IR surveys find fewer power-law AGNs as 
the smaller flux errors can result in larger x 2 values for 
the power-law fit. Thus, a source identified as a power- 
law AGN in a shallow survey may not be identified in a 
deeper survey as the more precise photometry are found 
to deviate from a power-law. Additionally the power-law 
sample is a very small subsample of the total X-ray pop- 
ulation, a major problem in selecting a large statistical 
population of AGNs. 

Donley et al. (2012) suggest adopting a 10% uncer- 
tainty floor on all flux measurements. Increasing the 
uncertainty to this level in our samples increases the 
number of selected Power-law objects, but the number 
density does not monotonically increase across our dif- 
ferent survey depths. While the power-law technique ap- 
pears to be a reliable method of selecting AGNs (a very 
high fraction of the power-law AGNs have X-ray coun- 
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Figure 3. X-ray area sensitivity curves that track the area which 
is sensitive for a given X-ray flux. We show the individual X-ray 
area sensitivity curves for the ~1 Ms (orange squares), ~100 ks 
(purples pluses), and ~10 ks (light blue crosses) samples. The 
vertical colored lines indicate our chosen X-ray flux limit for each 
survey depth, which is set at the flux which we recover 10% or 
more of the expected sources (see Table 1). The solid black line 
shows the All sample, which combines all of our X-ray fields that 
probe above the flux limit for each field. 

terparts, indicating they are true AGNs), it only selects 
a very small subsamplc of the AGN population (a tiny 
fraction of the X-ray sample and a much smaller sam- 
ple than identified with the Stern ct al. or Donley et al. 
selection techniques). Given the small sample size and 
the complex dependence on the depth of the IR data, we 
choose not to consider the power-law IR-AGN selection 
technique any further in the analysis of this paper. 

3.4. X-ray Selection 

We define samples, of X-ray selected AGNs based on 
detection in the hard (2—10 keV) energy band in our com- 
piled catalogs (see Section 2.3). A hard-band selected 
sample includes both X-ray unabsorbed and moderately 
absorbed AGNs, but misses Compton-thick sources. We 
do not include a soft-band selected sample even though 
it probes to fainter X-ray fluxes (providing a larger sam- 
ple) because it is more biased towards the detection of 
unabsorbed sources (with equivalent hydrogen column 
densities N H < 10 22 cm -2 ). 

Selecting in a single band also ensures we can accu- 
rately track (and correct for) the effects of variation in 
the X-ray sensitivity in our observations. Within an in- 
dividual Chandra or XMM-Newton pointing, variations 
due to the vignetting and point spread function (PSF) 
of the telescope lead to variations in the X-ray flux limit 
that is reached across the field of view. Furthermore, 
the exposure times between different pointings can vary 
within an X-ray survey, leading to substantial variations 
in sensitivity. The minor differences in the X-ray spec- 
tral shapes at 2—10 keV have a minimal effect on the 



sensitivity and the differences are mainly due to the X- 
ray telescope. Furthermore, at the X-ray fluxes we are 
probing (/ 2 -iok e v ~ > x 10~ 16 erg s~ f cm~ 2 ), the X- 
ray point source number counts are dominated by AGNs 
(Georgakakis et al. 2008; Lehmer et al. 2012). Thus, by 
comparing observed point source counts as a function of 
flux to the expected number based on the intrinsic Log N 

— Log S relation, we can estimate the fraction of our sur- 
vey area that is sensitive to X-ray AGNs of a given X-ray 
flux. 

In Figure 3 we show X-ray area curves (the area of 
a survey that is sensitive to sources above a given X- 
ray flux) for our various fields. We calculate these area 
curves using the ratio of the number of hard X-ray de- 
tected sources in the parent X-ray source catalog and the 
predicted number of sources based on the X-ray Log N 

— LogiS 1 relation of Georgakakis et al. (2008), using bins 
of 0.2 dcx width. We fix the area curve at the total area 
coverage above the lowest flux where >90% of the pre- 
dicted sources are detected to remove the effects of low 
source numbers at bright fluxes. To reduce the noise due 
to small numbers of sources at the faintest X-ray fluxes, 
we set a minimum flux threshold in each X-ray survey 
at the point where at least 10% of the expected number 
of sources are detected. These flux limits are shown in 
Figure 3 and given in Table 1. 

Similar to our IR observations, we group our X-ray 
fields into three different "X-ray survey depths" : 

1. — 1 Ms, which includes our deepest, 2 Ms X-ray 
data in the CDFS field. 

2. -100 ks, including the E-CDFS, COSMOS, the 
deeper regions of the XMM field. 

3. —10 ks, including ESI field and the remainders of 
the COSMOS and XMM fields. 

We give nominal hard band flux limits for each field and 
survey depth in Table 1. While there are slight differ- 
ences in the exact X-ray flux limit for each individual 
survey we show the surveys at one of the three "X-ray 
survey depths" for simplicity. However, we note that this 
flux limit is only reached over a small fraction (> 10%) of 
the total area covered by our fields with a given "X-ray 
survey depth" . 

Finally, we calculate an "X-ray weight" for each hard 
band X-ray detected source based on the ratio of the to- 
tal area coverage to the value of the area curve at the 
flux of the source. We use these X-ray weights to cor- 
rect observed number densities of X-ray sources to the 
intrinsic number density. A faint X-ray source will be 
given a large weight to account for the small area that is 
sensitive to faint fluxes. For the total X-ray sample, the 
X-ray weights ensure we recover the intrinsic X-ray Log N 

— Log S as measured by Georgakakis et al. (2008) that 
was originally used to determine the area curves. How- 
ever, the X-ray weights are also crucial when comparing 
between X-ray and IR-AGN samples (see Section 4). 

In Section 6, we further limit our X-ray detected sam- 
ple to objects with an observed Lx > 10 42 erg s _1 , calcu- 
lated from the hard band flux assuming a photon index 
of r = 1.9. This conservative limit ensures our sample 
contains objects where the X-ray flux is dominated by 
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emission due to an AGN, rather than star-formation pro- 
cesses in the galaxy. 

When comparing the X-ray AGN sample with the 
IRAC AGN samples we require that both samples lie 
within the overlapping window functions. Likewise for 
Section 6 we require that the objects fall within the re- 
gion with IRAC, X-ray and PRIMUS coverage. For these 
studies, we recompute our X-ray weights using the over- 
lapping areas only. 

4. SAMPLE COMPARISONS 

In this section we examine how the total numbers of 
sources and the overlap between our IR-AGN and X-ray 
AGN samples depend on the depths of the X-ray and 
IRAC data. In Section 4.1, we compare the IR-AGN 
and X-ray AGN selection techniques with respect to the 
bivariate IRAC and X-ray flux space. In Section 4.2, 
we determine the surface number densities of sources se- 
lected by the different techniques as a function of both 
IRAC and X-ray flux. In Section 4.3, we investigate how 
the overlap between the samples depends on the IR and 
X-ray survey depths. 

4.1. X-ray and IR Flux Comparison 

Our first step is to examine the distribution of sources 
in the bivariate space of observed IRAC 5.8 fim flux 
and hard-band X-ray flux. In the center panel of Fig- 
ure 4 we show contours (black) that trace the distribu- 
tion of all sources (across our differing survey depths) 
with both X-ray and IRAC detections. The contours 
contain 30%, 50%, and 80% of these sources; the remain- 
ing 20% of sources outside these contours are shown by 
black squares. We also show the distribution of sources 
that are also identified as IR-AGN using the Stern et 
al. (red) and Donley et al. (blue) selection techniques. 
Approximate IRAC and X-ray flux limits for the differ- 
ent survey depths are shown by the orange vertical and 
horizontal lines. We find that the majority of objects 
tend to be infrared bright when they are likewise X-ray 
bright. However, there is a population of objects that are 
infrared bright, but are X-ray faint, but very few that are 
X-ray bright but infrared faint. This is expected as that it 
suggests that it is easier to have a Compton-thick source 
that obscures the X-rays, but much harder to obscure the 
infrared photons. Many of these infrared-bright, X-ray 
faint sources are not selected by either the Stern et al. 
or Donley et al. IR-AGN selection techniques and thus 
may not be AGN or have dominate AGN components. 

The right panel of Figure 4 shows the number density 
of all IRAC and X-ray detected objects (black) , Stern et 
al. IR-AGN (red), and Donley et al. IR-AGN (blue) as 
a function of their observed hard band X-ray flux. Like- 
wise the bottom panel shows the number density of the 
same samples as a function of their IRAC 5.8 fxm flux. 
The number density plots in each of the sub panels are 
limited to include selected sources above the flux limit 
in each field and normalized by the area that is sensitive 
to that flux. The number densities of either Stern et al. 
or Donley et al. IR-AGN sample do not increase toward 
fainter fluxes as steeply as the total number density of 
objects detected. 



[5.8] Magnitude 
20 16 




flJBm Ny] 

Figure 4. The observed IRAC 5.8 /im and hard-band X-ray flux 
number-density distributions. In the center panel, we show the 
sample of objects detected with four bands of IRAC and a hard 
band X-ray detection. The Stern et al. IR-AGN selected sam- 
ple is shown with red squares, and the Donley et al. IR-AGN se- 
lected sample is shown with blue triangles. The approximate SWIRE, 
COSMOS, and GOODS depth IR surveys limits are shown with verti- 
cal orange lines, and approximate X-ray survey depths are shown 
with horizontal orange lines. The right side panel shows the total 
number of IRAC and X-ray detected of objects per logarithmic 
X-ray flux bin per area sensitive to that flux for the total (black), 
Stern et al. IR-AGN selected (red) and Donley et al. IR-AGN 
selected (blue) . The lower panel shows the total number of X-ray 
and IR-AGN detected objects per logarithmic IRAC flux bin per 
area sensitive to that flux for the respective samples. We include 
the X-ray completeness correction weights which corrects for the 
variation in X-ray flux limit in both side panels. 

4.2. X-ray and IR-AGN Surface Densities 

Our next goal is to determine how the number of 
sources identified as AGN using a given selection tech- 
nique varies as a function of IR flux. In the top panels 
of Figure 5 the solid black line shows the surface number 
density (number per square degree per logarithmic flux 
interval) of sources as a function of IRAC 5.8 /jm flux for 
the Stern et al. (left) and Donley et al. (right) IR-AGN 
samples. At the brightest fluxes (/5.8pm > 100 [iJy] we 
combine all of the different fields to estimate the sur- 
face number density. At fainter fluxes we only use the 
fields with sufficient IR depth (see Table 1 for flux lim- 
its). For both the Stern et al. and Donley et al. IR- 
AGN samples the number densities increase rapidly with 
decreasing flux with an approximately power-law form, 
although the slope of the power-law may flatten below 
/5.8pm ~ 100 fiJy for the Donley et al. sample. The 
dashed lines show the surface number densities of the 
Stern et al. or Donley et al. IR-AGN sample that is also 
detected in our X-ray data for our three different X-ray 
survey depths: ~10 ks (red), ~100 ks (green), and ~1 
Ms (blue). These surface number densities are corrected 
for the variation in X-ray flux limit over the field (the 
X-ray incompleteness) by applying the weights described 
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Figure 5. IR-AGN selected number density and fractions of combined IR-AGN samples as a function of IRAC fs.Siim flux for the Stern 
et al. (left column) and Donley et al. (right column) IR-AGN selection techniques. The top panels show the number density (black solid 
line) of the Stern et al. (left) and Donley et al. (right) IR-AGN selection method for the entire sample. The dashed lines show the number 
density for individual X-ray comparison depth samples: shallow ~10 ks (red), medium ~100 ks (green), and deep ~1 Ms (blue) samples are 
shown with Poisson error bars. For objects with measured X-ray fluxes we use the X-ray sensitive area derived from the X-ray completeness 
weights rather than the intersected IR and X-ray area. The total sample (solid black line) may dip below an individual X-ray survey depth 
sample due to the X-ray completeness corrections and small number statistics, but the black line is consistent with the individual X-ray 
survey depth lines. The bottom panels show the X-ray detected fraction of IR color selected AGN for each of the different X-ray survey 
depth bins. 



in Section 3.4. 

The bottom panels show the fraction of the IRAC sam- 
ples that are X-ray detected for each of the X-ray survey 
depths. These fractions vary significantly as a function 
of both IRAC flux and the depth of the X-ray data. At 
bright IR fluxes (/ 5 .8 M m > 100 fj,Jy) we find >70% of 
the IR-AGN samples are detected in X-ray data of ^100 
ks to ~1000 ks depth, whereas with the shallow X-ray 
data (~10 ks) this fraction is closer to 50%. At fainter 
IR fluxes the X-ray fraction decreases, although the fall 
off is highly dependent on the depth of the X-ray data. 
Figure 7 (left) directly compares the X-ray fractions of 
the Stern et al. and Donley et al. IR-AGN samples as 
a function of IR flux combining the different X-ray sur- 
vey depths. At bright IR fluxes ~60% of both IR-AGN 
samples are X-ray detected; at fainter IR fluxes a slightly 
higher fraction of the Donley et al. sources are identified 
in the X-ray data, although we note that the Donley et 
al. selection identifies a smaller sample of AGN. 

Our next step is to determine how the number of iden- 
tified IR-AGN and X-ray AGN sources depends on the 
X-ray flux. In the top panels of Figure 6 we show the sur- 
face number density of X-ray selected sources as a func- 
tion of the X-ray flux. We only include fields that probe 



to sufficient X-ray depth (see Table 1 and Figure 2) in our 
estimates of the surface number density at any given X- 
ray flux. We apply our X-ray weights to account for the 
variation in X-ray depth within a field (see Section 3.4); 
only a limited fraction of the combined X-ray and IRAC 
area may be sensitive to a given X-ray flux, even above 
our nominal X-ray flux limits. The X-ray surface number 
densities follow the double power-law distributions pre- 
sented by Georgakakis et al. (2008). The dashed lines 
show the surface number densities of X-ray sources that 
also satisfy the Stern et al. (left column) and Donley 
et al. (right column) IR-AGN selection criteria in fields 
reaching three different IR survey depths: SWIRE (red), 
COSMOS (green) and GOODS (blue). These number densi- 
ties are also corrected for variation in X-ray sensitivity 
within a field. The bottom panels show the fraction of 
the X-ray sample that satisfy the IR selection criteria for 
the three IR survey depths for each of the IR-AGN sam- 
ples. At bright X-ray fluxes (/x~10 -14 ergs _1 cm -2 ), 
a fairly high fraction of the X-ray sources (~60%) are 
identified as AGN using either the Stern et al. or Donley 
et al. selection criteria, although substantially increas- 
ing the depth of the IR data does not appear to increase 
this fraction. At fainter X-ray fluxes the fraction of X-ray 
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Figure 6. The number-density and fraction of the X-ray detected sample that is selected using a IR-AGN selection technique. The top 
panels show the number density of all X-ray detected objects (black solid line) and X-ray with IRAC detected objects(grey dotted linos) 
as a function of hard band X-ray. The dashed lines show the number density of the Stern et al. (left) and Donley et al. (right) IR-AGN 
samples for the shallow (SWIRE), medium (COSMOS), and deep (GOODS) depth IR surveys. Each of these dashed lines can be above the total 
number density of IR-AGN and X-ray AGN due to individual sample areas and small numbers at the lowest bins. We remove bins with 
objects less than 4 objects, which happen at the largest X-ray fluxes. The bottom panels show the IR-AGN fraction of X-ray AGN sample 
of the total number density for each individual IRAC depth survey. Objects with measured X-ray fluxes have been normalized by the 
sensitive area for that flux measurement rather than the intersected IR and X-ray area. 
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Figure 7. The X-ray detected fraction of IR-AGN (left panel) 
and the IR-AGN selected fraction of the X-ray AGN sample (right 
panel) for all field depths. The Stern et al. (solid red squares) 
and Donley et al. (solid blue triangles) IR-AGN selected fraction 
are for all samples which contains objects above the limits for that 
sample. Objects with measured X-ray fluxes have been normalized 
by the sensitive area for that flux measurement rather than the 
intersected IR and X-ray area. 



sources identified with the IR criteria reduces. Increasing 



the depth of the IR data from SWIRE to COSMOS recovers 
a higher fraction, but further increasing from COSMOS to 
GOODS depths does not appear to significantly increase 
the IR fraction at a given X-ray flux, especially for the 
Donley et al. selection criteria. The more restrictive 
Donley et al. selection criteria generally identifies a lower 
fraction of the X-ray sources as AGN at any given X-ray 
flux compared to the Stern et al. selection (see Figure 7, 
right panel, for a direct comparison of the fractions with 
COSMOS X-ray and IRAC depth data). 

With either of the IR-AGN selection techniques, the 
total number densities and the overlap between the IR- 
AGN and X-ray AGN samples depends on the survey 
depth. In the left panel of Figure 7 we show the X-ray 
detected fraction of the Stern et al. (red) and Donley et 
al. (blue) IR-AGN samples as a function of the IRAC 
flux. For sources fainter than /5.8pm ~ 100 fiJy IRAC 
flux, the Stern et al. fraction is smaller than the Don- 
ley et al. fraction, and very similar for sources above 
h.8nm ~ 100 fiJy. Below / 5 .8 Mm ~ 100 fiJy the Stern et 
al. selection identifies a larger number of objects as AGN 
than the Donley et al. selection (see Figure 8), which 
drives this fraction down. In the right panel of Figure 7, 
we show the Stern et al. (red) and Donley et al. (blue) 
IR-AGN fraction of the X-ray detected sample. The Don- 



10 



Mcndez et al. 



Icy ct al. IR-AGN selection identifies a lower fraction of 
X-ray detected objects for most of X-ray fluxes compared 
to the Stern et al. IR-AGN selection. 

4.3. Overlap of X-ray and IR-AGN samples 

In this section we examine in more detail the extent of 
the overlap between the IR-AGN and X-ray AGN sam- 
ples and how this depends on both the IR and X-ray 
survey depth. In Figure 8 (left panels) we show the frac- 
tion of the IR-AGN samples that are X-ray detected as a 
function of both the X-ray and IRAC flux limits for the 
different survey depths. The top-left panel shows the 
fraction of Stern et al. IR-AGN that are X-ray detected; 
the bottom-left panel shows the fraction of Donley ct al. 
IR-AGN that are X-ray detected. These fractions vary 
smoothly over the range of IR and X-ray flux limits be- 
tween 47% and 90% for both Stern et al. and Donley et 
al. IR-AGN samples. For comparable depth surveys in 
IR and X-ray (roughly the diagonal through these panels) 
70%-80% of the IR-AGNs are detected at X-ray wave- 
lengths. Increasing the depth of the IR data (moving 
to the left) reduces the fraction with X-ray detections; 
increasing the depth of the X-ray data (moving down) 
increases the fraction with X-ray detections. For our 
shallowest IR data we find a very high fraction (90%) of 
the IR-AGNs are detected using our deepest X-ray data. 
Generally the X-ray detected fraction of the Donley et 
al. IR-AGN samples is higher than for the Stern et al. 
IR-AGN samples, particularly for deeper surveys. 

In the two right-hand panels of Figure 8 we show the 
fraction of the X-ray AGN samples that are also identi- 
fied as Stern et al. IR-AGNs (top-right panel) or Donley 
et al. IR-AGNs (bottom-right panel). For both selection 
techniques, the fraction varies smoothly between ^4% 
and ^53%, being smallest for shallow IR and deep X-ray 
samples (bottom-right) and largest for deep IR and shal- 
low X-ray samples (top- left). The Donley et al. IR-AGN 
selection technique identifies a lower fraction of the X-ray 
sample than the Stern et al. IR-AGN selection technique. 
Comparing with Figure 6, we find that these low fractions 
are due to the large numbers of faint X-ray AGNs that 
are not selected by either IR-AGN technique, even with 
very deep IR data. 

In Figure 9 we use Venn diagrams to examine the sur- 
face densities and overlap of the samples selected using 
our three AGN selection techniques as a function of both 
X-ray and IR survey depth. The number density of the 
different sub-samples are represented by circles located at 
the IR and X-ray flux limits for the different surveys. The 
size of the circles are proportional to the surface number 
densities; the overlap between the samples are shown in 
Venn diagram form. The black cross indicates the X-ray 
and IR flux limits corresponding to the surveys used for 
each individual Venn diagram. The number density of 
X-ray AGN (green circles) increases rapidly as the X-ray 
depth increases (moving down in Figure 9). Likewise, 
for deeper IRAC surveys we find a larger number density 
of both Stern et al. and Donley et al. IR-AGN sam- 
ples. The overlap between the AGN samples varies sub- 
stantially depending on the IR and X-ray survey depth. 
X-ray selection identifies a large population of sources 
at all survey depths that are not identified with the IR- 
AGN selection techniques, even with extremely deep IR 
data. The Donley et al. IR-AGN selection technique 



generally identifies a subset of the Stern et al. IR-AGN 
sample for all depths. At our shallowest IR survey depth 
(SWIRE: right hand side of Figure 9) the IR-AGN samples 
overlap almost completely, while in deeper IR surveys an 
additional population is identified using the Stern et al. 
selection criteria. At the shallowest IR survey depth a 
high fraction of the IR-AGN are X-ray detected: up to 
90% with our deepest X-ray data. Deeper IR data identi- 
fies additional IR-AGNs, including populations that are 
not identified at X-ray wavelengths. 

5. CONTAMINATION OF IR-AGN SELECTION 

In this section we investigate whether either of the IR- 
AGN samples are subject to contamination from non- 
AGN galaxies. In Section 5.1, we examine the redshift 
distributions of the IR-AGN and X-ray AGN samples 
(using spectroscopic redshifts from PRIMUS at z < 1.2) 
and find clear evidence for contamination of the Stern et 
al. IR-AGN sample at specific redshifts. In Section 5.2, 
we use MIR templates to characterize the contamination 
from non-AGN host galaxies that fall into the Stern et al. 
IR-AGN selection criteria compared to the Donley et al. 
IR-AGN selection criteria and X-ray detected AGNs at 
z < 1.2. Finally in Section 5.3 we extend our template- 
based analysis to higher redshifts (z ~ 2—3), outside the 
range covered by the PRIMUS data, to investigate the 
contamination of the IR-AGN samples by star-forming 
galaxies at these redshifts. 

5.1. Investigating contamination with PRIMUS 
Redshifts 

In this section we compare the redshift distributions 
of our three AGN samples to investigate whether the se- 
lection techniques identifying the same underlying AGN 
population. In Figure 10, we show the number per square 
degree of Stern ct al. IR-AGN (red hatched), Donley 
et al. IR-AGN (blue hatched) and X-ray AGN (green 
solid) selected objects with PRIMUS redshifts. We re- 
strict our samples to objects with secure PRIMUS red- 
shifts between 0.2 < z < 1.2, where PRIMUS provides 
reliable redshifts for sources where the optical light is 
dominated by either the AGN or galaxy light and there 
is enough volume to construct a statistical sample. For 
the X-ray detected sample we ensure that the sample is 
not contaminated by X-ray detected star-forming galax- 
ies by restricting to sources with Lx > 10 42 erg s . We 
apply the weights calculated in Section 3.4 to correct for 
the varying depth of the X-ray data but we do not cor- 
rect for the fraction without PRIMUS redshifts. The 
fraction of objects without PRIMUS redshifts should be 
independent of selecting with either an IR-AGN or X- 
ray AGN technique. The Stern ct al. and Donley et 
al. distributions have been scaled up by a factor of 4 
to show their relative redshift distributions compared to 
the completeness-corrected X-ray sample. 

The Donley et al. IR-AGN sample and X-ray AGN 
sample number-density distributions peak around z ~ 
0.6, which matches the peak of the entire PRIMUS 
galaxy sample, although the Donley et al. distribution 
does not fall off as quickly for higher redshifts compared 
to the X-ray AGN sample. We find that the redshift 
distribution of the Stern et al. IR-AGN sample has an 
unexpected shape with two peaks, at z ~ 0.3 and z <~ 
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Figure 8. The X-ray detected fraction of IR-AGN selection techniques (left column) and the color selected fraction of X-ray AGN 
detections (right column) for each IR-AGN sample. The fractions are shown at the IRAC 5.8 /im and X-ray flux limits with background 
points shown with darker grey background for lower fractions, and lighter squares for higher fractions. The upper panels show the Stern et 
al. IR-AGN selected fractions, whereas the lower panels show the Donley et al. IR-AGN selected fractions. 



1.1, which differs significantly from that of the Donley 
et al. IR-AGN and X-ray AGN samples or the overall 
PRIMUS galaxy sample. To investigate this unusual dis- 
tribution, in Figure 11 we present IRAC color-color plots 
for the entire sample with IRAC detections in all four 
bands separated into four distinct redshift ranges. We 
over-plot the Donley et al. IR-AGN (blue triangle) and 
X-ray AGN (green diamond) selected samples. Gener- 
ally the Donley et al. IR-AGN sample is a subset of the 
Stern et al. IR-AGN sample, whereas the X-ray AGN 
sample includes objects both inside and outside of the 
Stern et al. color wedge. The Stern et al. IR-AGN sam- 
ple includes many other objects that are neither X-ray 
detected nor Donley et al. selected. We see that the ma- 
jority of these objects enter in at low redshifts (z ~ 0.3, 
left most panel) and high redshifts (z ~ 1.0, right most 
panel). These objects mostly enter into the wedge at 
low ([3.6] — [4.5]) color and appear to be consistent with 



scatter from the larger galaxy population (the dark red 
points in Figure 11). We thus conclude that the Stern 
et al. IR-AGN sample suffers from significant contam- 
ination by normal, non-AGN galaxies at redshifts z ~ 
0.3 and z ~ 1.1, resulting in the bimodal peaks in the 
redshift distribution seen in Figure 10. 

5.2. Comparison with IR SED Templates 

In this section we investigate the trends in the previous 
section using empirical galaxy and AGN SED templates. 
Following Donley et al. (2008), and Donley et al. (2012), 
we redshift the Polletta et al. (2008) IR templates and 
examine the changes to the observed IRAC colors. The 
Polletta et al. (2008) IR template set consists of four qui- 
escent galaxy templates, eight star-forming galaxy tem- 
plates, and six AGN templates, based on the observed 
SED of a representative sample of local galaxies and 
AGN. We convolve the redshiftcd SED with each IRAC 



12 



Mendez et al. 



10 



14 



X 
Z3 



o 10 

CO 



o 
a3 



15 



10" 16 h 



~1 1 — I — I — I — 

Stem 



~ i 1 — i — r 



Donley 

X-ray 

= 250 Number /deg 2 



XMM: 10-50ks + SWIRE 





ES1 :40ks + SWIRE 



COSMOS: 40ks + SCOSMOS 




CDFS: 
2Ms + SIMPLE 



CDFS: 2Ms + SWIRE 



CDFS: 
2Ms + GOODS 



10 1 



10 2 



/5.8,1m Ny] Flux Limit 



Figure 9. Venn diagrams of the number density and overlap between the Stern et al. (red), Donley et al. (blue), and hard band 
X-ray (green) selected AGN samples. The area and overlap of each circle is proportional to the surface number density of the sample and 
respective sample overlap. The different samples are located at the X-ray and IRAC flux limits shown as a black cross. Objects with X-ray 
detections have been normalized by the sensitive area for that flux, and included if they fall within the X-ray and IRAC limits for that 
sample. The X-ray and IRAC limit locations are plotted as a black plus at the center of the X-ray circle and Donley et al. selection circles, 
respectively. See Table 3 for individual fraction values. 



filter to predict the observed IRAC colors. In Figure 12, 
we show the predicted colors in four redshift bins for 
all the star-forming galaxy (blue dots), quiescent galaxy 
(red squares) and AGN (black crosses) templates. The 
lines track the change in color within the redshift bin for 
a given template. 

Generally star-forming templates change the most in 
this color-color space, due to the polycyclic aromatic hy- 
drocarbons (PAH) features in the SED. They start off 
with a very red ([5.8] — [8.0]) color that moves blue- ward 
due to the 8.0 /im PAH feature redshifting out of the 
8.0 \xm IRAC filter. The vertical ([3.6] -[4.5]) motion of 
star- forming templates is dominated by the 3.3 /im PAH 
feature moving from the 3.6 /xm to the 4.5 fim IRAC 
filter 14 . Together, both of these features can move the 

14 The normalization of the 3.3 jj,m feature in the templates 



star-forming templates into the Stern et al. wedge at 
both low redshift (z ~ 0.3) and high redshift (z ~ 1.1). 
The star-forming templates that come into the Stern et 
al. wedge are not selected by the Donley et al. IR-AGN 
criteria. 

For quiescent templates (red squares), we find a much 
simpler motion through the IRAC color-color space. 
Comprised of older stars, these templates arc gener- 
ally dominated by the Rayleigh- Jeans tail with a IRAC 
power-law slope of a ~ 2 located close to the origin. The 
SEDs arc generally featureless, and the motion is dom- 
inated by the 1.6 /-tm stellar bump which redshifts into 

may be too high relative to the continuum, which results in redder 
colors than are seen in the PRIMUS sample. However analyzing 
the detailed correspondence between star-forming galaxy templates 
and the IRAC colors of real galaxies is beyond the scope of this 
paper. 
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Figure 10. The redshift number-density distributions for the 
X-ray(grecn solid), Stern et al. (red hatched) and Donley et al. 
(blue hatched) selected samples for all objects with secure PRIMUS 
redshifts in all fields. The Stern et al. and Donley et al. IR- 
AGN distributions have been scaled up by a factor of four to easily 
compare their distributions to the X-ray sample. 

the 3.6 yon IRAC filter. These templates are generally 
excluded from either the Donley et al. or Stern et al. 
IR-AGN samples at these redshifts. 

The predicted AGN template colors also show very lit- 
tle motion in this space. Their featureless, red power-law 
(a <~ -0.5 - -3) SEDs tend to lie close to the power-law 
line. Both Donley et al. and Stern et al. IR-AGN selec- 
tion techniques pick up these AGN/Quasar (QSO) tem- 
plates easily at z < 1.2. 

The behaviors of the star-forming, quiescent, and AGN 
templates are consistent with the trends seen in the data 
in the previous section and explains the peaks in the 
Stern et al. IR-AGN redshift distribution. The Stern 
et al. IR-AGN sample is contaminated by star-forming 
galaxies at both low (z ~ 0.3) and high (z ~ 1.1) redshifts. 
The Donley et al. IR-AGN selection technique identifies 
AGN dominated templates at z < 1.2. 

The addition of a small AGN contribution (5%-10%) 
may move star-forming and quiescent SEDs into the 
Stern et al. wedge, but such sources are not selected by 
the Donley et al. criteria (e.g. Donley et al. 2012). Thus 
the Stern et al. selection may identify a more complete 
sample of AGN but it is hampered by a larger contami- 
nation by star-forming galaxies at particular redshifts. 

5.3. Extending Templates to Higher Redshifts 

In this section, we extend the galaxy and AGN tem- 
plates to higher redshifts (z > 1.2) beyond the range 
of PRIMUS to investigate the contamination of the IR- 
AGN samples. We track the templates out to z = 3 and 
compare to the full sample with IRAC detections thus 
including objects without spectroscopic redshifts. In Fig- 
ure 13 we show the density of objects in IRAC color-color 
space for the Donley et al. IR-AGN (top- left panel) , and 
X-ray AGN (top-right panel) samples, and for the entire 
sample (All) with IRAC detections (bottom-left panel). 
In the bottom-right panel we show the behavior of the 
templates from z = to z = 3. 

The Donley et al. IR-AGN sample (top-left panel) 



is centered around the power-law locus with the vast 
majority of the sample lying within the Stern et al. 
wedge. It spans the range of AGN/QSO templates (seen 
in the lower-right panel) over the entire redshift range 
(z = 0-3). 

In comparison the X-ray AGN sample (top-right panel) 
not only spans the range of the Donley et al. sample, 
but also extends beyond the Donley et al. IR-AGN sam- 
ple and the Stern et al. wedge. At ([3.5] -[4.5]) ~ 0.75 
and ([5.8]-[8.0]) > 0.5 (left of the Stern et al. wedge) 
there is a spur of X-ray sources that do not have a se- 
cure PRIMUS redshift and are not selected with the Don- 
ley et al. IR-AGN method. Comparing this spur with 
the template track panel, this region of objects appears 
to match both high redshift (z > 2.5) quiescent galaxy 
(square) and star- forming galaxies (cross). It is unclear 
whether these objects contain an X-ray bright AGN and 
arc missed by the Donley et al. and Stern et al. methods, 
or are merely X-ray-detected high-redshift star-forming 
galaxy contaminants. If we assume that these objects 
have a z = 3 redshift, the hard band X-ray luminosity 
would be distributed around L\ > 10 44 ' 6 erg s _1 , which 
is suggestive that they are highly luminous AGN and not 
high-redshift star-forming galaxy contaminants. The ob- 
jects below the Stern et al. wedge generally match qui- 
escent galaxy templates with z < 1.5 and star-forming 
templates with 0.7 < z < 1.5. This population has been 
shown in Section 5.1 to be low redshift galaxies with 
PRIMUS redshifts (z ~ 0.7) containing genuine AGN 
(Lx > 10 42 erg s _1 ) that are not identified by cither IR- 
AGN selection technique. 

The small population of X-ray detected sources with 
very red colors (([5.8]-[8.0]) > 2) to the bottom right 
of the Stern et al. wedge are generally identified as 
low-redshift (z < 0.4) star-forming galaxies with low X- 
ray luminosities (Lx < 10 42 erg s _1 ) and thus the X-ray 
emission may be attributed to star-formation rather than 
AGN activity. This population is easily identified and 
not selected by either the Stern et al. or Donley et al. 
IR-AGN selection techniques. 

For the entire sample of IRAC detected objects (All, 
bottom left panel), the bulk of the objects lie below 
the Stern et al. wedge, concentrating at colors of 
([3.6]-[4.5]) < 0.0 and ([5.8]-[8.0]) > 0.1. The All sam- 
ple contains four-bands IRAC detected objects from all 
of our fields. The colors of these objects correlate with 
predicted star-forming and quiescent galaxy templates 
colors at 0.7 < z < 1.5 There is also a large density of 
sources to the left side of the Stern et al. wedge, like the 
z > 2.5 spur of X-ray detected objects discussed earlier. 
The colors of these objects match both quiescent and 
star-forming galaxies SED templates at high-redshift (z 
~ 2.5). A low fraction (~5%) of these objects are high- 
redshift high-luminosity X-ray detected AGN, but the 
majority of these objects are not X-ray detected. These 
objects are not identified by the Donley et al. selection 
technique, but extend into the Stern et al. wedge con- 
taminating it with high-redshift galaxies. 

Some of these sources in the Stern et al. wedge could 
have a low luminosity AGN, but the AGN signature is 
being masked in the MIR by the dominant galaxy light. 
Such AGN could be identified using other techniques (e.g. 
X-ray, or via their optical spectra), but the true AGNs 
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Figure 11. The IRAC color-color space for sources with PRIMUS redshifts in four redshift bins between 0.2 < z < 1.2. In each panel 
the Stern et al. wedge is outlined in orange with the power-law —0.5 < a < —3 line also shown in orange within the wedge. The Donley 
et al. IR-AGN sample (blue triangles) and X-ray AGN sample (green diamonds) is shown above all other objects (dark red circles) in each 
redshift bin. We find that the bimodal redshift peaks of the Stern et al. IR-AGN sample are due to a large number of objects entering the 
bottom of the wedge at both z ~ 0.3 and z ~ 1.1. 
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Figure 12. The IRAC color-color spaces for galaxies and AGN templates in four redshift bins between 0.2 < z < 1.20 to match the 
PRIMUS redshift sample. In each panel we show the redshifted colors of templates in that redshift bin (lines). The Stern et al. wedge 
is outlined in orange with the power-law —0.5 < a < —3.0 line also shown in orange within the wedge. Quiescent galaxies(red square), 
star- forming galaxies (blue cross) and QSO/AGN (black diamond) redshift template tracks begin in each redshift bin at the filled symbol, 
and continue for the extent of that redshift bin. Generally we find the Donley et al. and Stern et al. IR-AGN techniques select QSO/AGN 
templates. The Stern et al. IR-AGN technique additionally selects star-forming templates at both z ~ 0.3 and z ~ 1.1. Quiescent galaxy 
templates are not selected by either at these redshifts. 



are probably a minority of the contaminating population. 
The dominant galaxy light suggest that they are not gen- 
erally broad-line QSOs, but could host narrow-line AGN. 

In summary, the Stern et al. IR-AGN technique selects 
a large population of AGN but is contaminated by low-, 
intermediate-, and high-redshift galaxies (z ~ 0.3, z ~ 
1.1, and z > 2.5). The Donley et al. IR-AGN selection 
technique generally is free of these galaxy contaminants, 
selecting just the AGN population. The X-ray detected 
sample has colors that extend beyond the Stern et al. 
and Donley et al. selected samples, identifying a popu- 
lation of AGN that are not detected by cither IR-AGN 
selection technique and are not star-forming galaxy con- 



taminants. 

6. AGN AND HOST GALAXY REST-FRAME 
PROPERTIES 

In this section we investigate the AGN and host galaxy 
rest-frame properties of our IR-AGN and X-ray AGN 
samples. In Section 6.1, we create subsamples with 
PRIMUS spectroscopic redshifts in order to study rest- 
frame properties. In Section 6.2, we measure the AGN lu- 
minosity and hardness ratio distributions for the IR and 
X-ray AGN samples. We then investigate how the X-ray 
luminosities vary across the IRAC color-color space. In 
Section 6.3, we compare the host galaxy colors and stellar 
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Figure 13. Comparison of the Donley et al. IR-AGN and X- 
ray AGN samples to the entire IR-AGN sample and the IR SED 
templates in IRAC color-color space. In the Donley ct al. (top 
left), X-ray (top right) and All (bottom left) panels we indicate 
the density of IRAC detected objects by the greyscale color. We 
use a linear greyscale where black corresponds to the maximum 
density of objects in a given panel. The Stern ct al. wedge is 
outlined in orange with the power-law —0.5 < a < —3.0 line also 
shown in orange within the wedge. In the bottom right panel, 
we show the iRAC colors of the templates from z = to z = 3 
(indicated by the color coding) for star-forming galaxies(cross at 
z = 0), quiescent galaxies (square at z = 0), and AGN (diamond 
at z = 0). The templates indicate that the Stern et al. IR-AGN 
sample is contaminated by galaxies at high redshifts(z>2), which 
is confirmed by our observational dataset (the over-density in the 
All panel to the left of the Stern et al. wedge). 

masses for the IR and X-ray AGN samples to investigate 
differences in the host galaxy properties. 

6.1. The PRIMUS Redshift Sample 

To investigate AGN and host galaxy restframc proper- 
ties we require a PRIMUS spectroscopic redshift. This 
reduces our AGN sample sizes and imposes an optical 
apparent magnitude limit. To be targeted by PRIMUS 
an object must have an optical magnitude brighter than 
i ~ 23. The fraction of sources in our AGN samples that 
are brighter than this optical limit varies from —30% to 
—70% depending on the depth of the corresponding IR or 
X-ray data. Of these sources that are brighter than the 
PRIMUS targeting limit, -20% to -60% were targeted 
in PRIMUS (not all objects to the apparent magnitude 
limit can be targeted due to slit collisions). The fraction 
of targeted AGN that have a high quality redshift varies 
from —75% to —95% depending on the depth of the IR 
or X-ray data, where the fraction is higher for shallower 
fields. This is comparable to the overall PRIMUS red- 
shift success rate. Some targeted sources for which we 
fail to measure a redshift lie at z > 1.2 where the B aimer 
break falls beyond the wavelength coverage of PRIMUS. 
Our redshift success range in the 0.2 < z < 1.0 range 
is likely to be higher. The fraction of our IR and X-ray 
selected AGN with PRIMUS rcdshifts is -20% to -70% 
depending on depth of the data. We find that the sub- 
sample with PRIMUS redshifts spans the full range of 
IR and X-ray fluxes in our IR and X-ray AGN samples. 



In the rest of this section we use these subsamples to 
compare the restframe properties of the IR and X-ray 
selected AGN to an optical limit of i — 23. 

6.2. Rest-frame AGN Properties 

In this section we investigate rest-frame properties re- 
lated to the AGN themselves. In Figure 14 we compare 
the IR and X-ray luminosities for the AGN samples. For 
the IR luminosity we use £3.6^1, defined to be vL v at a 
rest-frame wavelength of 3.6 fim, calculated in two man- 
ners: 1) interpolating the log of the fluxes using the two 
nearest IRAC channels, and 2) using template fit syn- 
thesized luminosities from K-correct (Blanton & Rowcis 
2007). We find no significant differences in the luminos- 
ity distributions between these two methods and use the 
former for the remainder of the paper. The solid green 
circles indicate PRIMUS AGN that have IR and X-ray 
detections and are above Lx > 10 42 erg s _1 . X-ray lumi- 
nosities are calculated from the 2 — 10 keV fluxes, us- 
ing an unobscured power-law with a photon index of 
r = 1.9. The light green arrows indicate upper limits on 
the 3.6 [im IR luminosity for X-ray detected AGN that 
lack IR detections. The IR luminosity upper limits are 
derived from the IR depth of the survey and the redshift 
of the object. The X-ray luminosity upper limits for the 
Stern et al. or Donley ct al. IR-AGN that are not X-ray 
detected are derived from the X-ray flux corresponding 
to 50% of the expected number of X-ray sources are de- 
tected in that field. For a more conservative 90% limit, 
the upper limits increase by roughly a factor of 2—4. Ob- 
jects that are selected by either Stern et al. or Donley 
et al. IR-AGN techniques are indicated by red squares 
or blue triangles, respectively. Generally, sources that 
are detected in both IR and X-ray have correlated lumi- 
nosities which roughly follows the Richards et al. (2006) 
mean QSO SED shown as an solid orange line. The corre- 
lation does contain significant scatter with more sources 
scattering to much lower X-ray luminosities. The sources 
that scatter to lower X-ray luminosities are found two or- 
ders of magnitude below the mean QSO SED 

We find that both IR-AGN color selection techniques 
tend to select sources with higher IR and X-ray lumi- 
nosities compared to the entire X-ray selected AGN sam- 
ple. A weighted KS test of the distributions of the X- 
ray luminosities, as shown in the top-right subpanel of 
Figure 14, rejects the null hypothesis that either the 
Stern et al. or Donley et al. IR-AGN are drawn from 
the same distribution as the X-ray selected AGN at the 
> 99.9% confidence level. X-ray selection identifies a 
large population of AGN with low luminosities (below 
Lx < 10 43 ' 5 erg s _1 ) that are not identified as cither 
Stern et al. or Donley et al. IR-AGN. 

These trends are also shown in Figure 15(a), which 
plots the subset of X-ray AGN detected in all four IRAC 
bands in IRAC color-color space, where the color of the 
points scale with X-ray luminosity. Low X-ray luminos- 
ity (L x < 10 43 - 5 erg s _1 ) AGN are scattered throughout 
the IRAC color-color space, whereas high X-ray luminos- 
ity (Lx > 10 43 ' 5 erg s _1 )AGN cluster near the power- 
law line within the Stern et al. wedge. For high X-ray 
luminosity (L x > 10 43 5 ergs" 1 ) AGN with IRAC de- 
tections, -46%±6% to -79%±4% lie within the spaces 
defined by the Donley et al. and Stern et al. IR-AGN 
selections, respectively. 
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Figure 14. X-ray vs. IR Luminosity comparison plot (top left) for the objects that are X-ray or IRAC detected. The Stern et al. selected 
objects are shown with a red square, the Donley et al. selected objects are shown with smaller blue square, and the X-ray detected objects 
with Lx > 10 42 erg s -1 shown in green filled squares. We show the X-ray detected objects with redshifts but not detected in the IRAC 
with light-green upper limits at the depth of the IR survey. We show the Stern et al. or Donley et al. IR-AGNs with redshifts but not 
X-ray detections with light red or blue upper limits at the X-ray flux corresponding to 50% of the expected number of X-ray sources are 
detected for each field. The IR and X-ray sample luminosities are correlated; the more luminous IR sources are also the more luminous 
X-ray sources. This correlation matches the Richards et al. (2006) mean quasar SED, which we show as an orange solid line. In the top 
right panel we show the X-ray luminosity number-density distribution for the sources in the previous panel. We show the Stern et al. 
distribution with a red hatched histogram, the Donley et al. distribution with a blue diagonal-hatched histogram, and the X-ray and (not) 
IR detected sample distribution in (light) green filled (arrows) circles. We do not show Stern et al. or Donley et al. IR-AGN that are not 
X-ray detected in this histogram. In the bottom left panel we show the residuals between the source luminosity and the Richards et al. 
(2006) mean quasar SED line for the X-ray and IR detected objects. For the distribution panels, objects with X-ray detections have been 
normalized by the sensitive area for that flux. The points arc shown with the same colors and styles as the top left panel. We find that 
the residuals generally scatter below the mean SED line with more sources having a higher IR luminosity than X-ray luminosity. 



In the bottom panel of Figure 14 we plot the residu- 
als X-ray to IR luminosity ratio from the Richards et al. 
(2006) mean Quasar SED shown in the upper panel (solid 
orange line). In this panel we have removed the AGN 
not detected in the IR. We find that there is a tail of 
sources with negative residuals such that their X-ray to 
IR luminosity ratio is below the Richards et al. (2006) 
mean quasar SED (solid orange) line. There are be- 
tween 3%±1% and 6%±3% of Stern et al. and Donley et 
al. sources with residuals less than —2.0; these sources 
are particularly X-ray under luminous. These sources are 
likely obscured AGN where the observed X-ray luminos- 



ity is suppressed due to absorption in the 2 — 10 keV 
X-ray band. 

We also examine the obscuration of the X-ray AGN us- 
ing hardness ratios (HRs). We convert the hard- and soft- 
band fluxes to equivalent on-axis Chandra count rates 
and calculate the HRs as (H-S)/(H+S) where H and S 
are the 2—10 keV and 0.5—2 keV count rates, respectively. 
Large HRs values (HR > 0.2) indicate significant absorp- 
tion in the hard X-ray band (Njj > 3x 10 21 cm~ 2 ) at z 
~ 0.6 (Hasinger 2008). In Figure 15(b) we show the HRs 
for the subset of X-ray AGN detected in all four IRAC 
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Figure 15. (a) IRAC color-color plot of the X-ray detected population where we scale the color of each object to the X-ray luminosity. 
We show the border of the Stern et al. wedge as a solid black line, and the power-law locus as a dashed black line. We find a greater than 
average number of high X-ray luminosity sources to be within the Stern et al. wedge, and more than average numbers of sources with low 
X-ray luminosity to be outside of the wedge. There is a population of medium X-ray luminosity(Lx ~ 10 43 ' 5 erg s — 1 ) sources that are 
not selected by the Stern et al. or Donley et al. techniques, (b) IRAC color-color plot of the X-ray detected population with object color 
scaled to the X-ray flux hardness ratio. Generally we find that Hard sources (HR~ 1) are not preferentially found within the Stern et al. 
wedge compared to the full sample. 
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Figure 16. Hardness ratio comparison for the Stern et al. (hori- 
zontal hatched red shading), Donley et al. (diagonal hatched blue 
shading), and X-ray (green filled shading) identified samples. The 
two panels show the X-ray hardness ratio split up into X-ray lumi- 
nous sources (Lx > 10 43 ergs -1 ) (bottom left), and X-ray faint 
sources (Lx < 10 43 ergs -1 ) (bottom right). In the bottom two 
panels, we find that the population of less luminous AGN sources 
(bottom left) that are missed by the IRAGN sources have higher 
hardness ratios compared to the X-ray sample and thus are more 
obscured. 
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Figure 17. IR luminosity (L^.e^m) number density distributions 
for sample intersections. In both panels we show the X-ray distri- 
bution in green. In the left panel we show the number density of 
Stern et al. sources that are X-ray detected (red solid) and the 
number density of Stern et al. sources that are not X-ray detected 
(red hatched). In the right panel we show the number density of 
Donley et al. sources that are X-ray detected (blue solid) and 
the number density of Donley et al. sources that are not X-ray 
detected (blue hatched). We do not include the X-ray complete- 
ness correction weights for cither IR selection. We find that the 
IR-AGN samples that are not X-ray detected have similar median 
luminosities compared to the X-ray detected sample. 
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bands in IRAC color-color space and find that lower HR 
(unobscured) AGN tend to fall within the Stern et al. 
wedge. This is consistent with IR-AGN samples identi- 
fying primarily more luminous AGN, which tend to be 
unobscured. Obscured AGN lie throughout the IRAC 
color-color space, with 48%±3% in the Stern et al. wedge 
and the remainder in the region where the galaxy light 
dominates the IR colors. 

We compare the HRs of the IR-AGN and X-ray 
AGN samples in a two X-ray luminosity bins to limit 
the large number of luminous IR-AGN sources from 
dominating the comparison. In Figure 16(a) we show 
the distributions of HRs for X-ray luminous sources 
(L x > 10 43 ergs" 1 ). The Stern et al. IR-AGN, Don- 
ley et al. IR-AGN, and X-ray detected AGN samples 
all span a wide range of HRs. The fraction of obscured 
sources (HR > 0.2) in the X-ray AGN sample (28%±7%) 
is slightly larger than the fraction in the Stern et al. 
(18%±10%) or Donley et al. (20%±17%), although the 
differences are not statistically significant. 

In Figure 16(b) we show the distributions of HRs for 
lower luminosity X-ray sources (Lx < 10 43 erg s _1 ) and 
find that the fraction of obscured sources (HR > 0.2) 
for the X-ray AGN sample (28%±5%) is slightly smaller 
than the fraction in the Stern et al. (36%±15%) or Don- 
ley et al. (67%±27%) AGN samples. Unlike at high 
X-ray (Lx > 10 43 erg s _1 ) luminosities, there is a hint 
at low X-ray (Lx < 10 43 erg s _1 ) luminosities that the 
IR-AGN samples identify a more obscured population; 
however, there are not statistically significant differences 
between the HR distributions in either luminosity range. 
Larger deep-IR samples are needed to confirm if this 
small population of obscured sources statistically differs 
from the X-ray sample. 

Finally, in Figure 17 we show the 3.6 fim IR luminos- 
ity number density distribution for the IR-AGN samples 
that are divided into X-ray detected (solid histograms) 
and not X-ray detected (hatched histograms). We also 
show a subset in green of the X-ray AGN sample that is 
detected in all IRAC four bands. For both the Stern et al. 
and Donley et al. AGN samples, we do not find any signif- 
icant differences in the Lir distributions for the sources 
with X-ray detections to those without X-ray detections. 
These latter sources could be missed in X-ray due to 
heavy obscuration, but nevertheless they have a similar 
distribution of IR luminosities and do not outnumber the 
subsamplc of the X-ray AGN sample detected in IRAC. 
The low number density of IR-AGN sources without X- 
ray detections may partially be due to having relatively 
deeper X-ray survey data compared to the IRAC survey 
data. Using the 3.6 /zm and X-ray luminosity from the 
Richards et al. (2006) mean quasar SED, we find that the 
depth of the X-ray data in our fields is typically 2-3 times 
deeper than the IRAC depth, for light that is dominated 
by luminous AGNs. Thus while IR-AGN samples are ef- 
ficient at identifying luminous AGN, it remains unclear 
if they identify a large population of obscured sources 
relative to X-ray AGN samples due to the differences in 
survey depths. 

6.3. Host Galaxy Properties 

In this section we investigate the optical rest-frame col- 
ors and stellar masses of the host galaxies for our IR- 



AGN and X-ray AGN samples with PRIMUS redshifts. 
For broad-line AGN the optical light is dominated by the 
AGN rather than the host galaxy light; this prevents an 
accurate estimate of the host galaxy properties. There- 
fore in this section we exclude sources classified as broad- 
line AGN, where the PRIMUS spectrum is better fit with 
\ 2 of at least 50 by a broad-line AGN template com- 
pared to any galaxy template. These broad-line AGN are 
a sizable fraction of the X-ray AGN (12%±1%), Stern 
et al. IR-AGN (16%±2%) and Donley et al. IR-AGN 
(32%±3%) samples. 

In the top panels of Figure 18, we plot the rest-frame 
optical (u — g) colors and stellar masses for the X-ray 
AGN (left), Stern et al. IR-AGN (center), and Donley et 
al. IR-AGN (right column) samples. We show the rele- 
vant AGN sample with black squares and PRIMUS galax- 
ies with gray contours (30%, 50%, and 80% contours) and 
gray points for individual galaxies outside of the 80% con- 
tour. In the left panel we show the X-ray AGN that are 
selected by either of the Stern et al. IR-AGN or Donley 
et al. IR-AGN selection techniques with green diamonds. 
In the center and right panels we show IR-AGN that 
are also X-ray detected in green diamonds. We use the 
stellar mass estimates of Moustakas et al. (2013), which 
are based on SED fits to multi-wavelength photometry 
from the UV to the near-IR. In the bottom panels of Fig- 
ure 18, we calculate the fraction of all PRIMUS galaxies 
that are detected in one of the AGN selection methods as 
a function of stellar mass for red and blue sources. We 
define each source as blue or red using the magnitude 
and redshift dependent cut from Aird et al. (2012): 

C = (u - g) - (0.671 - 0.031M g - 0.065z) (10) 

where M g is the absolute g-band magnitude and z is the 
redshift. To calculate the number density and account 
for the variation in survey depth, we use the X-ray com- 
pleteness corrections (X-ray sensitive area for a given X- 
ray flux) for the sources that are X-ray detected and the 
IR-Optical overlap area for all other sources. 

We find that the fraction of galaxies with detected X- 
ray AGN increases with stellar mass. The fraction of blue 
galaxies with X-ray AGN is a factor ^2 higher than for 
red galaxies (across all stellar masses). These trends are 
in agreement with Aird et al. (2012), where the increase 
with stellar mass is attributed to a selection effect. AGN 
in massive galaxies appear more luminous (for the same 
specific accretion rate scaled relative to the host stellar 
mass) and so are easier to detect. AGN selected by either 
Stern et al. IR-AGN or Donley et al. IR-AGN techniques 
follow the same basic trend, although both techniques 
generally identify a lower fraction of AGN in both red 
and blue galaxies at all stellar masses. The rate that 
the fraction of red galaxies with Stern et al. or Donley 
et al. AGN increases with stellar mass is much slower 
than for red galaxies with X-ray AGN. This lower rate of 
change for the IR-AGN selections suggests that they are 
less efficient at identifying AGNs in red, massive galaxies 
(M* > lO 10 ' 5 Mq). The fraction of galaxies with Stern 
et al. AGN is higher than the X-ray fraction at low 
masses which may be due to contamination of Stern et 
al. IR-AGN by star-forming galaxies. 

In Figure 19 (a) we estimate the specific accretion rate 
(A) of the X-ray detected sample as a function of the host 
stellar mass. Foiling Aird et al. (2012), we define the 
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Figure 18. Optical color versus stellar mass diagram for the X-ray AGN (left column), Stern et al. IR-AGN (center column) and Donley 
et al. IR-AGN (right column) samples. In the top panel of each column we show the rest-frame (u — g) optical colors on the vertical axis 
and the stellar mass estimates for PRIMUS galaxies (grey points and contours) and the AGN sample (black squares). In each panel we 
show the sources that are both X-ray detected and IR-AGN selected with green diamonds. The background colored hatching shows the 
approximate red and blue galaxies definition at z = and m g ~ 20. In the bottom panel of each column we show the fraction of detected 
AGN compared to all PRIMUS galaxies at a given stellar mass for red and blue hosts. The fraction is calculated in six stellar mass bins 
and is shown at the median AGN mass in that bin. 



specific accretion rate from the bolometric luminosity de- 
rived from the X-ray luminosity using the Hopkins et al. 
(2007) quasar bolometric corrections 15 . For the hard X- 
ray selected sample this is dominated by the emission 
from the AGN with possibly only minor contribution 
from the host galaxy. From the bolometric luminosity 
we calculate the specific accretion rate: 

, _ Lboi / 1 1 s 

^ Edd 1.3 x 10 38 erg s" 1 x 0.0024^ 

where is the host galaxy stellar mass and L^i is the 
bolometric luminosity. The X-ray sources with IRAC de- 
tections are solid green circles and the sources without 
IRAC detections are open green circles. Sources with 
IRAC detections that are selected by either Stern et al. 
or Donley et al. IR-AGN selection technique are outlined 
with red squares or blue diamonds. The vast majority of 
Donley et al. IR-AGN also also satisfy the Stern et al. 
IR-AGN selection criteria. This figure shows two popula- 
tions of X-ray sources that are not identified as IR-AGNs: 
those that are not detected by IRAC (see Section 4.3) 
and those that are detected by IRAC but do not meet 
either the Stern et al. or Donley et al. criteria. The 
former populations consists primarily of sources at low 

15 https://www.cfa.harvard.edu/ phopkins/Site/qlf.html 



stellar masses and high specific accretion rates. These 
sources are primarily found in the fields that have shal- 
lower IRAC data compared to their X-ray data (ESI, XMM, 
C0SM0S shallow ). 

The majority of the X-ray AGNs that are not Stern et 
al. or Donley et al. selected (but do have IRAC detec- 
tions: green solid symbols in Figure 19 (a)) have large 
stellar masses (A4* > 3x 10 10 Mq), low-to-moderate 
X-ray luminosities (Lx < 10 43 ergs _1 ), and low spe- 
cific accretion rates (A < 1%). Additionally, many 
of these sources have a moderate-to-large MIR to X- 
ray luminosity ratios (£3.6/™/ Lx> !)• In fact, their 
L3.Qp.ra/Lx ratios are larger than typically found for the 
IR-AGN population and for optically selected quasars 
(e.g. Richards et al. 2006). This suggests that stellar 
light from their massive host galaxies is dominating the 
emission in the IRAC bands. Indeed, for a stellar mass- 
to-light ratio of order unity (e.g. Bell & de Jong 2001) 
and L 3 . Sfim / Lx for a pure AGN SED (Richards et al. 
2006), a galaxy with a specific accretion rate of A = 1% 
would have comparable emission at 3.6 \xm from stellar 
light and from the central AGN. At lower specific ac- 
cretion rates, the host galaxy would then dominate the 
emission in the IRAC bands. This provides a natural 
explanation for why many X-ray AGN at lower specific 
accretion rates are not selected as IR-AGNs. 

In Figure 19 (b) we estimate the specific accretion rate 
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Figure 19. (a) Specific accretion rate (A) derived from the X-ray luminosity versus stellar mass. For sources that are detected in both 
IRAC and X-ray (filled dark green circles) we also show Stern et al. or Donley et al. IR-AGN sources as red squares or blue diamonds, 
respectively. The grey dashed line shows the approximate location of an Lx ~ 10 42 erg s — 1 source. We find many X-ray detected sources 
with high specific accretion rates (A > 1CT 1 ) are not IRAC detected, especially at low mass, (b) Specific accretion rate (A) derived from 
the L3.6 M m luminosity versus stellar mass. X-ray detected sources are filled green circles, with filled red squares or blue diamonds for Stern 
et al. and Donley et al. IR-AGN, respectively. Stern et al. and Donley et al. sources that are not X-ray detected are shown with open light 
red squares and light blue diamonds, respectively. There is a floor near 10~ 2 < A < 10~ 3 where galaxy light begins to dominate the MIR 
emission, (c) Specific accretion rate (A) distribution for the X-ray, Stern et al. and Donley et al. AGN samples for Ai t > 3x 10 10 A4q. 
The Stern et al. and Donley et al. IR-AGN techniques tend to identify higher specific accretion rate sources relative to the X-ray AGN 
sample. We include the X-ray completeness corrections for the X-ray detected sources. 



based on MIR luminosity. The specific accretion rates for 
these sources are less well constrained due to the possible 
addition of galaxy light to the observed MIR flux. To en- 
sure that the bolomctric luminosity distribution are sim- 
ilar for the X-ray and IRAC overlapping sample, we first 
use the Richards et al. (2006) mean quasar SED to esti- 
mate the Lx luminosity from the L^,Q^ m luminosity (see 
the orange line in Figure 14). We apply the Hopkins et al. 
(2007) hard X-ray bolomctric correction to the equivalent 
X-ray luminosity to calculate the specific accretion rate 
for these sources. We do not use the Hopkins et al. (2007) 
bolomctric corrections for the IR luminosity due to the 
overprediction of the bolometric luminosity for IR-AGN 
sources that are also X-ray detected. This method for es- 
timating specific accretion rates illustrates the floor near 



10~ 2 < A < 10~ 3 , where galaxy light dominates the MIR 
luminosity. 

In Figure 19(c) we show the specific accretion rate 
distribution for Stern et al. IR-AGN, and Donley et 
al. IR-AGN and X-ray AGN. We restrict this com- 
parison to AGN in host galaxies with stellar mass 
M* > 3x 10 10 Mq, due to incompleteness at lower stel- 
lar masses. For all AGN with an X-ray detection, we 
use the specific accretion rate derived from the X-ray lu- 
minosity. For the IR-AGNs that arc not detected in the 
X-rays we use their estimated specific accretion rate from 
their MIR luminosity. Some of the IR-AGNs with X-ray 
detections have significant absorption in the X-rays (e.g. 
Figure 14). For IR-AGNs with L 3 . efim /L x > 10, we use 
the specific accretion rates derived from their MIR lumi- 
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nositics. Wc note that this docs not have a significant 
effect on the lambda distributions in Figure 19(c). Both 
the Stern et al. and Donley et al. IR-AGN samples have 
specific accretion rate distributions that do not match 
the X-ray distribution at the > 99.7% confidence level 
from using a two-sided weighted KS test. In particu- 
lar, the X-ray sample includes a significant number of 
sources at low to moderate specific accretion rates that 
are not identified as IR-AGNs. These low- to moderate- 
specific accretion rate sources are preferentially found in 
red, high-mass galaxies, which explains the larger frac- 
tion of red host galaxies identified with X-ray AGN com- 
pared to IR-AGN selection techniques. 

7. DISCUSSION 

In this paper we compare AGN samples selected in 
the MIR to those selected using hard-band X-rays, in or- 
der to quantify both the overlap and uniqueness of each 
selection method as well as understand the AGN pop- 
ulations identified. We use the Stern et al. and Don- 
ley et al. IR-AGN selection methods (described in Sec- 
tion 3) to identify IR-AGN samples using Spitzer/IRAC 
data in four fields that have X-ray coverage from Chan- 
dra or XMM-Newton: CDFS, COSMOS, Elais-Sl (ESI), 
and the XMM-LSS field. We create IR-AGN and X-ray 
AGN samples to various depths in both the IR and X-ray 
data and find that the number density of AGN recovered 
with each selection method varies strongly as a function 
of depth, as does the overlap between the IR-AGN and 
X-ray AGN samples. We use PRIMUS spectroscopic red- 
shifts in these fields to study the AGN luminosities, host 
galaxy colors, and host stellar masses of the IR-AGN and 
X-ray AGN samples and investigate potential contamina- 
tion in IR-AGN selection. 

There are several advantages to our approach, com- 
pared to previous studies. We use large AGN samples (a 
total of 9,005 Stern ct al. IR-AGN, 3,569 Donley ct al. 
IR-AGN and 4,886 X-ray AGN) spanning multiple fields, 
which allows us to perform statistical comparisons as a 
function of joint parameters, while minimizing the effects 
of cosmic variance. We take advantage of the large num- 
ber of spectroscopic redshifts in these fields provided by 
PRIMUS (a total of 1,541 AGNs with PRIMUS redshifts 
between 0.2< z <1.2), so that we do not have to rely on 
photometric redshifts which have larger errors. We fur- 
ther study the overlap and uniqueness of IR-selected ver- 
sus X-ray selected AGN as a function of both the depth 
of the IR data and the depth of the X-ray data. With 
a large dynamic range in survey depth, we can test how 
the completeness, purity and contamination of the AGN 
selection depends on depth. Such tests are required to 
reconcile the often vastly different and seemingly conflict- 
ing results in the literature regarding IR and X-ray AGN 
selection. 

Another key difference in our methodology is that we 
use hard-band (2 — 10 keV) X-ray selection, as opposed to 
soft- or full-band, and we use weights for X-ray detections 
to account for incompleteness in the X-ray data. Hard- 
band X-ray selection is sensitive to both unabsorbed and 
moderately-absorbed AGN, though it will miss heavily 
absorbed, Compton-thick AGN. Soft (0.5 - 2 kcV) or 
full (0.5 — 10 kcV) band selection, by comparison, is bi- 
ased towards unabsorbed sources (with column densities 
of Nh < 10 22 cm -2 ). We account for the large positional 



uncertainties in the X-ray data by using the likelihood ra- 
tio matching technique to securely match the X-ray sam- 
ple to both the IR and optical samples. Finally, we ac- 
count for the varying sensitivity of the X-ray data across 
individual fields by applying statistical completeness cor- 
rections to all X-ray sources. This last correction is cru- 
cial to accurately quantify the number density of X-ray 
AGN as a function of depth and robustly compare with 
IR-AGN samples. We emphasize that this correction is 
not generally applied in the literature, where studies sim- 
ply compare the fraction of IR-AGN that are observed in 
X-ray data, without taking into account incompleteness 
in the X-ray data. 

7.1. How AGN Selection Varies with IR and X-ray 

Depth 

Wc first consider the fraction of IR-AGN that are de- 
tected in X-rays, as a function of both the IR and X-ray 
data depth. Generally, a higher fraction of Donley et 
al. IR-AGN are X-ray detected than in the Stern et al. 
IR-AGN sample, consistent with our findings above that 
the Stern et al. IR-AGN selection can suffer from con- 
tamination. At the relatively shallow IR depth of the 
SWIREsurvey, 60% of IR-AGN are X-ray detected when 
the X-ray depth is comparable to the IR depth. Using 
the deepest IR and X-ray data, the fraction of IR-AGN 
that are X-ray detected rises to ^70% to ~80%. How- 
ever, this fraction rises to ^90% when combining shallow 
IR and deep X-ray data, indicating that the vast major- 
ity of IR-AGNs are indeed X-ray emitting AGNs. 

As discussed above, IR-AGN selection methods iden- 
tify additional AGN not detected in X-rays, when the IR 
and X-ray data are of comparable depth. At shallow IR 
depths the Stern et al. and Donley et al. techniques re- 
cover very similar number densities of AGN with a large 
overlap between the populations. As IR depth increases, 
both the Stern et al. and Donley et al. selections iden- 
tify a larger population of AGNs. However, the number 
density recovered by the Stern et al. selection increases 
at a greater rate - in our deepest IR data Stern ct al. re- 
covers approximately twice the number density of AGN 
as the Donley et al. IR-AGN selection. As discussed in 
Section 5, much of the relative increase is likely due to 
contamination. For any depth IR and X-ray data, the 
addition of Stern et al. IR-AGN sample to X-ray AGN 
sample increases the total AGN sample by 7%— 27%, de- 
pending on the depth of the IR and X-ray data. In 
comparison, addition of Donley et al. IR-AGN selection 
to X-ray AGN selection increases the total AGN sample 
by 4% -20%. Both IR-AGN selection methods increase 
the total AGN samples more at shallower IR and X-ray 
depths. At all IR and X-ray depths studied here, X-ray 
AGN selection identifies a higher number density of AGN 
than either IR-AGN selection method. 

It is also extremely important to account for the well- 
known and easily characterized incompleteness of X-ray 
selection. The variation in X-ray sensitivity over a field 
can substantially reduce the number of IR-AGN with ob- 
served X-ray detections. It is vital to account for this 
effect to understand any true, underlying differences be- 
tween AGN populations 

At all IR and X-ray depths studied here we find 
a much higher fraction of IR-AGN arc detected in 
X-rays, compared to values in the literature. This 
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is due to the X-ray completeness corrections that we 
make, to account for the varying X-ray sensitivity 
within a field. If we do not apply these complete- 
ness weights, we find similar fractions to those in the 
literature. For example, in the XBootes survey with 
shallow IR (/5.8/jm ~ 51 fJ-Jy) and shallow X-ray 
(/x~10 -15 erg s^cm- 2 ) data, Hickox et al. (2009) find 
that 38% of their Stern et al. IR-AGN are X-ray de- 
tected. Using similar depth data if we do not apply X- 
ray completeness corrections we find that 37% of Stern 
et al. IR-AGN arc X-ray detected; however this fraction 
rises to 63% after applying the completeness corrections. 
Donley et al. (2007) use deep IR (/5.8pm ~ 14.5 [iJy) 
and medium-depth X-ray (/x~10 -15 erg s _1 cm~ 2 ) data 
in the CDF-North field to find that 33% of the Stern et al. 
IR-AGN are X-ray detected. Using similar depth data if 
we do not apply X-ray completeness corrections we find 
a similar fraction of 35% of Stern et al. IR-AGN are 
X-ray detected; however this fraction rises to 57% after 
applying the completeness corrections. Park et al. (2010) 
use very deep IR data (/5.8/xm ~ 6.3 fiJy) and medium- 
depth X-ray (/ x ~ 10~ 16 erg s^cm" 2 ) data in the EGS 
field and find that a comparatively low fraction of Stern 
et al. IR-AGN are X-ray detected. As their X-ray data is 
not as deep as their IR data, it is not surprising that they 
find a low fraction. We do not have comparably deep IR 
data, however at a similar X-ray depth using the deepest 
IR data we have, we find that at most 35% of IR-AGN 
should be X-ray detected, if X-ray completeness is not 
accounted for. This upper limit is consistent with their 
value. 

Thus, we find that the wide range of X-ray detected 
IR-AGN fractions reported in the literature can be ac- 
counted for by the varying depths of both the IR and 
X-ray data used in these studies. We also find that ac- 
counting for X-ray incompleteness, which is generally not 
done in the literature, increases this fraction by a fac- 
tor ^2 which results in 25% more of the IR-AGNs be- 
ing counted (statistically) as X-ray detected. This illus- 
trates how differing depths of data can lead to very differ- 
ent conclusions in prior studies of IR and X-ray AGNs. 
It is also extremely important to account for the well- 
known and easily characterized incompleteness of X-ray 
selection. The variation in X-ray sensitivity over a field 
can substantially reduce the number of IR-AGN with ob- 
served X-ray detections. It is vital to account for this 
effect to understand any true, underlying differences be- 
tween AGN populations. 

7.2. Contamination and Bias of IR-AGN Selection 

Much, if not most, of the difference between the Stern 
et al. IR-AGN sample and Donley et al. IR-AGN sample 
is due to contamination by galaxies that do not host an 
AGN (or at least, are not dominated by a luminous AGN 
in the MIR). Using PRIMUS redshifts, we study the evo- 
lution of the distribution of Stern et al. and Donley et al. 
IR-AGN and X-ray detected AGN in IRAC color-color 
space from z~0.2toz~1.2 and compare our results 
to the predicted evolution of star-forming and quiescent 
galaxy and AGN SED templates in this space. We find 
that the Stern et al. IR-AGN sample is contaminated at 
z ~ 0.3 by star-forming galaxies and at z ~ 1.1 by qui- 
escent galaxies. This is reflected in the fraction of Stern 



et al. IR-AGN that are not detected in X-rays; for the 
full Stern et al. IR-AGN sample with PRIMUS redshifts, 
spanning 0.2 < z < 1.2, this fraction is 27%±4%, while 
at 0.2 < z < 0.45 it is 76%±10% and at 1.0 < z < 1.2 it 
is 86%±13% when using the X-ray completeness weights. 

Extending our analysis of the galaxy and AGN tem- 
plates to higher redshift, beyond the reach of PRIMUS, 
we find that the Stern et al. IR-AGN samples will be 
contaminated by high-redshift star-forming galaxies. We 
do not find any galaxy contamination in the Donley et 
al. IR-AGN sample; this is reflected by the lack of evo- 
lution in the fraction of Donley et al. IR-AGN that 
are not X-ray detected. Using the Donley et al. IR- 
AGN sample as a baseline, we find that galaxy contam- 
ination in the Stern et al. IR-AGN selection is signifi- 
cant only for IR surveys deeper than the SWIRE limit 
of /5.8pm ~ 100 [iJy. 

Overall, we find that the Donley et al. IR-AGN selec- 
tion is less complete than the Stern et al. IR-AGN se- 
lection. This is demonstrated by the fact that a smaller 
fraction of X-ray AGN are also Donley et al. IR-AGN. 
Additionally there are correctly identified AGNs within 
the Stern et al. IR-AGN sample that are not found by 
the Donley et al. IR-AGN selection. However, the Don- 
ley et al. IR-AGN selection is more reliable than the 
Stern et al. IR-AGN selection, in that a higher fraction 
of Donley et al. IR-AGN are also X-ray detected, and 
the Donley et al. IR-AGN selection does not suffer from 
galaxy contamination. The Stern et al. and Donley et al. 
IR-AGN samples demonstrates the importance of finding 
a balance between minimizing both contamination and 
incompleteness . 

Both the Donley et al. and Stern et al. IR-AGN selec- 
tions appear to be biased, however, in that they select 
high luminosity AGN. This is seen in the X-ray lumi- 
nosity distributions of the IR-AGN samples compared 
to the X-ray AGN sample, where the IR-AGN selec- 
tion techniques preferentially sample the population of 
Lx > lO 43 - 5 ^ erg s _1 AGN. These high luminosity sources 
have IRAC colors similar to power-law AGN and are eas- 
ily detected with shallow IR surveys. 

7.3. Uniqueness of IR-AGN Selection 

Even with a relatively deep X-ray survey and shallow 
IR data, we find that there is at least ~10% of the IR- 
AGN population which is not detected in X-rays even 
after correcting for the variable X-ray sensitivity. For 
samples of comparable X-ray and IR depth, this fraction 
is typically ^20% to ^30%. This reflects the uniqueness 
of IR-AGN selection, in that these are sources that are 
only identified in the IR and would typically be missed 
in X-ray surveys. What kind of AGN does IR selection 
uniquely identify? These AGN could in theory be missed 
by X-ray selection either because they are intrinsically 
less X-ray luminous, due to a lower accretion rate onto 
the SMBH, or they could be obscured by a high column 
density of gas and dust, limiting the ability of even deep 
X-ray surveys to detect them. For the population of IR- 
AGN that are X-ray detected, as discussed above we find 
that IR-AGN selections techniques preferentially select a 
population of luminous AGN with Lx > 10 4 erg s" 1 . 
However, there is a small fraction (~8%) of IR-AGN with 
very high IR to X-ray luminosity ratios, implying heavy 
obscuration. Therefore up to ~20% of IR-AGN are likely 
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sources with moderate to heavy obscuration. 

In comparing the hardness ratio distributions of IR and 
X-ray selected AGN, we find no statistically significantly 
difference in the distributions for X-ray AGN with higher 
or lower X-ray luminosity. This comparison does not 
account for the ~10% of IR-AGN sources that are not 
X-ray detected. Thus while the majority of IR-AGN are 
also X-ray AGN, the IR selection techniques adds a small 
but important population of obscured sources that are 
missed even with the deepest X-ray surveys. 

7.4. Host Galaxies 

With the exception of broad-line sources, the optical 
SEDs of the IR-AGN and X-ray AGN samples are dom- 
inated by light from the host galaxy, which allows us to 
determine galaxy properties for these sources. Generally 
the Stern et al. and Donley et al. IR-AGN samples fol- 
low the same trends as the X-ray population, in that IR- 
AGN and X-ray AGNs do not populate vastly different 
host galaxy populations. The galaxy hosts span a similar 
range of stellar masses (10 10 < M*/ M Q < 10 115 ) and 
optical colors compared to the PRIMUS parent galaxy 
sample. The Donley et al. IR-AGN and X-ray AGN 
samples identify similar fractions of AGN in blue, star- 
forming hosts and red, quiescent hosts. These fractions 
of red and blue AGN hosts are similar to the red and 
blue fractions for all PRIMUS galaxies. 

Hickox et al. (2009) find that X-ray AGN are preferen- 
tially found in galaxy hosts in the green valley, the mini- 
mum of the optical color bi-modality between the red se- 
quence and the blue cloud (Martin et al. 2007), whereas 
Stern et al. IR-AGN are typically in blue host galaxies. 
We do not agree with these results, which did not correct 
for X-ray incompleteness and contamination in the Stern 
et al. sample. Most recent papers (e.g. Silverman et al. 
2009; Xue et al. 2010; Aird et al. 2012) find that the X- 
ray AGN fraction is highly stellar mass dependent. Sim- 
ilar to Nandra et al. (2007), Silverman et al. (2009) and 
Aird et al. (2012) we find that at a given stellar mass the 
fraction of galaxies with X-ray AGN is higher in the blue 
cloud. In particular for the green valley Nandra et al. 
(2007) using the ^200 ks X-ray data in the Extended 
Groth Strip to find that 13% of X-ray detected AGN 
have green host galaxies, which is consistent with our 
X-ray AGN sample (9%) using a similar green valley def- 
inition. If X-ray completeness correction weights are not 
applied, the fraction of green host galaxies with X-ray 
AGN rises only slightly, from 9% to 11%. 

For the Stern et al. IR-AGN we generally find equal 
fractions of sources with blue or red host galaxies. There 
is a slightly larger fraction with blue hosts when com- 
pared using stellar-mass matched samples, but this is 
mainly due to blue star-forming host galaxies contaminat- 
ing the sample rather than differences in the underlying 
population. The Donley et al. IR-AGN sample is more 
evenly distributed between red and blue host galaxies at 
fixed stellar mass. 

We do find that the fraction of galaxies hosting an 
AGN increases with galaxy stellar mass. The fraction of 
massive (M* > 10 11 Mq) galaxy hosts is larger for the 
X-ray AGN sample (78%±2%) compared to the Stern 
et al. IR-AGN (52%±2%) or Donley et al. IR-AGN 
sample (69%±4%). Conversely the fraction of less mas- 
sive (10 9 < M*/ M Q < 10 10 ) galaxy hosts for the X- 



ray AGN sample (7%±1%) is smaller either the Stern et 
al. IR-AGN sample (21%±2%) or the Donley et al. IR- 
AGN sample (11%±3%). However, these differences are 
minor compared to the similarity of the host galaxy stel- 
lar mass distribution between any of the AGN selection 
techniques. 

From the host galaxy stellar masses and the estimated 
AGN bolometric luminosities we infer the specific accre- 
tion rates, which gives us a valuable measure of the 
intrinsic differences of the selected AGN populations. 
While IR-AGN samples are generally found in similar 
host galaxies as X-ray AGN, we find that IR-AGN have 
relatively high specific accretion rates. The X-ray AGN 
sample includes a large number of sources at low to mod- 
erate specific accretion rates that are not identified as 
IR-AGNs. These sources reside in red, high-mass galax- 
ies that have large IR to X-ray luminosities suggesting 
that they are not being identified by the IR-AGN tech- 
niques due to a bright galaxy component in the MIR 
disguising the low accretion rate AGN. They contribute 
to the large rise in the fraction of red massive galaxies 
that have a X-ray AGN. This accounts in part for the sig- 
nificant population of X-ray AGN that is not identified 
by either IR-AGN technique. 

8. CONCLUSIONS 

In this paper we compare the completeness, contami- 
nation, overlap and uniqueness of Spitzer /IRAC and X- 
ray identified AGN. We quantify the differences in X- 
ray AGN and IR-AGN selection techniques due to IR 
and X-ray survey depths, X-ray obscuration, and sur- 
vey completeness to assess the usefulness of these tech- 
niques to identify obscured AGNs. We use Spitzer /IRAC 
data, XMM-Newton and Chandra X-ray data at multiple 
depths to construct the largest sample of IR and X-ray se- 
lected AGN to date. We focus on the Stern et al. (2005) 
and Donley et al. (2012) selection techniques which se- 
lect AGNs that dominate the MIR with a red power-law 
SED. We statistically compare the IR-AGN populations 
against X-ray selected AGN samples. The combination 
of multiple X-ray and IRAC depth surveys in four fields 
(CDFS, ESI, COSMOS, and XMM) allows us to study 
the overlap and uniqueness of the AGN selections as a 
function of both IR and X-ray depth; having multiple 
fields also reduces cosmic variance. This gives us a large 
sample with 9,005 Stern et al. IR-AGN, 3,569 Donley et 
al. IR-AGN and 4,886 X-ray AGN. Characterizing the 
variation of the samples as a function of survey depth 
also allows us to probe the bias and contamination of 
individual identification techniques. 

We take advantage of the ~1,500 secure PRIMUS red- 
shifts in these fields to probe the intrinsic properties of 
the AGNs, including luminosity, specific accretion rate, 
and hardness ratio, as well as properties of the host galax- 
ies, including color and stellar mass. We compare the 
AGN properties of the IR-AGN versus X-ray AGN sam- 
ples and compare their host galaxies to the full PRIMUS 
galaxy population. The main results from our work are 
as follows: 

1. IR-AGN selection identifies predominantly lumi- 
nous AGNs, with L^.Q^m ~ 10 44 5 ergs -1 and 
L x ~ 10 43 - 5 ergs" 1 . X-ray AGN selection iden- 
tifies a larger population of lower luminosity and 
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lower accretion rate AGNs, where the host galaxy 
dominates the MIR emission. These AGNs are not 
identified using IR-AGN selection. 

2. IR-AGN selection does not identify a substantial 
population of obscured AGNs relative to X-ray 
samples. However, ~10% of IR-AGN are not de- 
tected in extremely deep X-ray data and could be 
very heavily obscured. 

3. Stern ct al. IR-AGN selection is contaminated 
by non-AGNs at specific redshifts (z ~ 0.3, z ~ 
1.1, and z > 2.5). The level of contamination de- 
pends on depth and is not significant at shallow 
(e.g. SWIRE) IR depths. Donley et al. IR-AGN 
selection is not contaminated. 

4. IR-AGN and X-ray AGN samples both preferen- 
tially identify AGNs in massive galaxies. They fur- 
ther both identify AGN in red and blue host galax- 
ies, with a similar ratio of red to blue galaxies. The 
host galaxy stellar masses and colors are therefore 
quite comparable between IR and X-ray AGN se- 
lection. 

5. Both the Stern et al. and Donley et al. IR-AGN 
samples identify AGNs with high specific accretion 
rates relative to the X-ray detected AGN sample. 
The low and moderate accretion rate AGN iden- 
tified only in X-rays have large stellar mass host 
galaxies, which dominate the MIR SED. 

We find that the majority of current X-ray and IR sur- 
veys in cosmological fields are often not well matched in 
depth, in that the X-ray data is generally significantly 
deeper than the IR data. For example, a pure AGN 
SED shows that GOODS-depth IRAC data and ~20 ks 
X-ray data will probe roughly the same intrinsic AGN 
luminosity. Therefore, while existing deep X-ray surveys 
ensure that we have a fairly complete X-ray AGN sam- 
ple, deeper IR surveys are needed to build a comparably 
complete IR AGN sample. Moreover, larger shallow and 
wide X-ray surveys can be used to accurately compare to 
the luminous IR-AGN samples from these shallower and 
wider IR surveys. 

IR and X-ray AGN selections appear to identify fairly 
similar AGN populations in similar high stellar mass 
galaxies, with no strong preference for either red or blue 
host galaxies. There is a large overlap between these 
selections, though IR selection preferentially identifies 
brighter AGN. However IR-AGN selection techniques do 
identify a small population (<~10%) that is not identified 
in extremely deep X-ray surveys and could be very heav- 
ily obscured. A combined IR and X-ray AGN selection 
will identify a more complete sample than cither selection 
alone, including both heavily obscured AGNs and AGNs 
that arc under-luminous relative to their host galaxies. 
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Table 1 

IRAC and X-ray Flux Limits. 



Field X-ray Depth IR Depth X f AC /r 5 - 8 /T t- -.r -f -2! 

J Limit [fiJyl Limit [erg s 'cm z \ 



CDFS 


~2Ms 


GOODS 


20.0 


l.BxlO -18 






COSMOS 


44.7 


1.8xl0 -18 






SWIRE 


100.0 


1.8xl0 - 18 




~200ks 


GOODS 


20.0 


7.1xl0~ 16 






COSMOS 


44.7 


7.1X10" 16 






SWIRE 


100.0 


7.1xl0~ 16 


COSMOS 


~160ks 


COSMOS 


44.7 


7.1xl0- ift 




~40ks 




44.7 


2.8xl0~ 15 


ESI 


~40ks 


SWIRE 


100.0 


2.8xl0- ib 


XMM 


~50-100ks 


SWIRE 


100.0 


7.1xl0" it) 




~10-50ks 




100.0 


5.6xl0~ 15 



Table 2 

Field Areas and Overlap. 



Area [deg 2 ] Overlap Area [deg 2 ] 



Field 


X-ray Depth 


IR Depth 


X-ray 


IRAC 


PRIMUS 


IRAC 
X-ray 


PRIMUS 
IRAC 


PRIMUS 

X-ray 


PRIMUS 
IRAC 
X-ray 


CDFS 


~2Ms 


GOODS 


0.12 


0.06 


2.90 


0.06 


0.06 


0.12 


0.06 






COSMOS 


0.12 


0.44 


2.90 


0.12 


0.21 


0.12 


0.12 






SWIRE 


0.12 


7.24 


2.90 


0.12 


2.57 


0.12 


0.12 




~200ks 


GOODS 


0.31 


0.06 


2.90 


0.06 


0.06 


0.18 


0.06 






COSMOS 


0.31 


0.44 


2.90 


0.30 


0.21 


0.18 


0.18 






SWIRE 


0.31 


7.24 


2.90 


0.31 


2.57 


0.18 


0.18 


COSM09 


~160ks 


COSMOS 


1.00 


2.61 


1.18 


1.00 


1.18 


0.81 


0.81 




~40ks 




2.16 


2.61 


1.18 


2.16 


1.18 


1.18 


1.18 


ESI 


~40ks 


SWIRE 


0.64 


6.03 


1.03 


0.64 


1.03 


0.59 


0.59 


XMM 


~50-100ks 


SWIRE 


1.07 


8.28 


3.41 


1.03 


3.37 


0.91 


0.88 




~10-50ks 




5.43 


8.28 


3.41 


4.07 


3.37 


1.30 


1.30 



Table 3 

Raw Numbers and Number Densities of IR-AGN and X-ray AGN samples in each field with overlapping IRAC and X-ray Data. 



Field 



Limits Weighted Number / deg 2 (Raw Number) 

X-ray Depth IR Depth X-ray a Stern Donley 



Overlap Weighted Number / deg 2 (Raw Number) 

, Stern 
Donley „ , 

Donley 



Stern 
Donley 



Stern 
X-ray a 



X-ray a 



X-ray a 



CDFS 


~2Ms 


GOODS 


5055 (176) 


2477 (115) 


1118 (53) 


1118 (53) 


1769 (74) 


928 (42) 


928 (42) 






COSMOS 


4912 (245) 


1135 (91) 


709 (59) 


709 (59) 


996 (74) 


627 (49) 


627 (49) 






SWIRE 


4912 (245) 


336 (33) 


269 (26) 


269 (26) 


303 (29) 


244 (23) 


244 (23) 




~200ks 


GOODS 


2204 (81) 


1598 (84) 


729 (38) 


729 (38) 


752 (35) 


436 (21) 


436 (21) 






COSMOS 


2077 (418) 


961 (240) 


521 (139) 


521 (139) 


787 (187) 


436 (113) 


436 (113) 






SWIRE 


2073 (421) 


349 (96) 


244 (69) 


244 (69) 


317 (86) 


218 (61) 


218 (61) 


COSMOS 


~160ks 


COSMOS 


2067 (1176) 


972 (766) 


556 (443) 


549 (436) 


717 (510) 


438 (324) 


433 (320) 




~40ks 




750 (1104) 


700 (1279) 


408 (744) 


403 (735) 


399 (628) 


256 (416) 


253 (411) 


ESI 


~40ks 


SWIRE 


740 (162) 


345 (133) 


259 (99) 


257 (98) 


250 (72) 


182 (50) 


180 (49) 


XMM 


~50-100ks 


SWIRE 


2167 (1083) 


348 (306) 


244 (226) 


242 (224) 


295 (251) 


198 (178) 


198 (178) 




~10-50ks 




371 (892) 


272 (902) 


214 (717) 


212 (709) 


151 (409) 


111 (296) 


111 (296) 



For distributions where we require X-ray detections, we use the X-ray completeness weights. 
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Table 4 

Overlap Fraction of Stern IR-AGN and X-ray AGN. 



CD ^ 

as o p, a) v 

g -g cn ^-co oco 

^ >> h :^ « 3 R ft += ft 

xl &t g§ la 

13 S « ^ 



o " o 



a g£ o 



S 13 feo 



Field X-ray Depth IR Depth P-, £ fc£ PU T3 P -8 



CDFS 


~2Ms 


GOODS 


71.4%±0.9% 


35.0% ±0.7% 


30.7%±0.6% 


12.3%±0.4% 


57.0%±0.7% 






COSMOS 


87.7%±1.0% 


20.3%±0.6% 


19.7%±0.6% 


2.8%±0.2% 


77.5%±0.6% 






SWIRE 


90%±2% 


6.2%±0.3% 


6.1%±0.3% 


0.7%±0.1% 


93.2%±0.4% 




~200ks 


GOODS 


47%±1% 


34%±1% 


24.7%±0.8% 


27.7%±0.8% 


47.6%±0.9% 






COSMOS 


82%±1% 


38%±1% 


35%±1% 


7.7%±0.6% 


57%±1% 






SWIRE 


91%±2% 


15.3%±0.8% 


15.0%±0.8% 


1.5%±0.3% 


83.4%±0.8% 


COSMOS 


~160ks 


COSMOS 


74%±1% 


35%±1% 


30.9%±1.0% 


11.0%±0.6% 


58%±1% 




~40ks 




57%±2% 


53%±2% 


38%±1% 


29%±1% 


33%±1% 


ESI 


~40ks 


SWIRE 


72%±2% 


34%±2% 


30%±2% 


11%±1% 


59%±2% 


XMM 


~50-100ks 


SWIRE 


85%±2% 


13.6%±0.7% 


13.3%±0.7% 


2.4%±0.3% 


84.3%±0.8% 




~10-50ks 




55%±3% 


41%±3% 


31%±2% 


25%±2% 


45%±2% 



Table 5 

Overlap Fraction of Donley IR-AGN and X-ray AGN. 



CO H 171 o CO 



< 73 ^ "-5 

^ J< «e 
-gB £cg 



r. 



g3 2o bo 

S£ £H 

Field X-ray Depth IR Depth PU S :g p_ tj X 8 



CDFS 


~2Ms 


GOODS 


83%±1% 


18.4%±0.5% 


17.7%±0.5% 


3.6%±0.3% 


78.7%±0.6% 






COSMOS 


88%±1% 


12.8%±0.5% 


12.6%±0.5% 


1.6%±0.2% 


85.8%±0.5% 






SWIRE 


91%±2% 


5.0%±0.3% 


4.9%±0.3% 


0.5%±0.1% 


94.6%±0.3% 




~200ks 


GOODS 


60%±2% 


19.8%±0.8% 


17.5%±0.8% 


11.7%±0.6% 


70.8%±0.9% 






COSMOS 


84%±2% 


21.0%±0.9% 


20.2%±0.9% 


3.9%±0.4% 


75.9%±0.9% 






SWIRE 


89%±2% 


10.5%±0.7% 


10.4%±0.7% 


1.2%±0.2% 


88.4%±0.7% 


COSMOS 


~160ks 


COSMOS 


79%±2% 


21.2%±0.9% 


20.0%±0.9% 


5.4%±0.5% 


74.5%±0.9% 




~40ks 




63%±2% 


34%±2% 


28%±2% 


17%±1% 


55%±2% 


ESI 


~40ks 


SWIRE 


70%±3% 


25%±2% 


22%±1% 


9%±1% 


68%±2% 


XMM 


~50-100ks 


SWIRE 


81%±3% 


9.1%±0.6% 


8.9%±0.6% 


2.1%±0.3% 


89.0%±0.7% 




~10-50ks 




52%±3% 


30%±2% 


23%±2% 


22%±2% 


55%±2% 



